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DYNAMIC RESPONSES OF RAILROAD CAR MODELS 


TO VERTICAL AND LATERAL RAIL INPUTS 

By John L. Sewall, Russell V. Parrish, 
and Barbara J. Durling 
Langley Research Center 

SUMMARY 

Simplified dynamic models were applied in a study of vibration in a high-speed 
railroad car. The mathematical models used were a four-degree-of -freedom model 
for vertical responses to vertical rail inputs and a ten-degree-of -freedom model for 
lateral responses to lateral or rolling (cross-level) inputs from the rails. Elastic 
properties of the passenger car body were represented by bending and torsion of a uni- 
form beam. Rail-to-car (truck) suspensions were modeled as spring-mass-dashpot 
oscillators. Lateral spring nonlinearities approximating certain complicated truck 
mechanisms were introduced. The models were excited by displacement and, in some 
cases, velocity inputs from the rails by both deterministic (including sinusoidal) and 
random input functions. Results were obtained both in the frequency and time domains. 

Solutions in the time domain for the lateral model were obtained for a wide variety 
of transient and random inputs generated on-line by an analog computer. Variations in 
one of the damping properties of the lateral car suspension gave large fluctuations in 
response over a range of car speeds for a given input. This damping coefficient was 
significant in reducing lateral car responses that were higher for nonlinear springs than 
for linear springs for three different inputs. Results showed that knowledge of car stiff- 
ness properties is needed in order to locate system resonances which need to be avoided. 

INTRODUCTION 

The recent introduction of high-speed rail passenger service between metropolitan 
areas has stimulated considerable interest in the effects of vibration on passenger com- 
fort. Vibrations are caused by track irregularities, and reduction of these vibrations at 
high speeds is a particularly difficult task. 

References 1 to 5 report some of the analytical research directed toward under- 
standing the parameters influencing the generation and transmission of vibrations in rail 
systems. In all these studies, significant motions and structural characteristics in highly 


complicated suspension mechanisms are represented by simple spring-mass-damper 
models. In reference 1, accelerations and forces transmitted from the rails were studied 
for ranges of spring constants and damping coefficients of a three-degree-of -freedom 
model of vertical motion that included railbed flexibility. Lateral accelerations and 
forces were studied for a seven-degree-of -freedom model that was also used to explore 
yaw stability. Reference 2 is a study of the responses of a simplified model of an auto 
ferry transport to transient inputs. In reference 3, a four- degree- of- freedom model is 
utilized to show the advantages of a new suspension design over a conventional design in 
reducing accelerations transmitted to the car. Reference 4 is a vibration-minimization 
study of a three-degree-of -freedom model of an electrically powered railroad car subject 
to vertical sinusoidal and random inputs from the rails. In reference 5 a four-degree- 
of -freedom vertical model and a ten-degree-of -freedom lateral model were used to 
examine a wide variety of parameter variations, including the effects of car bending 
stiffness and spring nonlinearity on car response. 

Topics not sufficiently considered in previous studies are lateral responses to both 
transient and random inputs, optimization that includes the effects of truck-suspension 
(equalizer) damping and stiffness, and damping effects in models which include spring 
nonlinearities. None of these studies (including the present one) allows for rigid-body 
pitch motion of the car; consequently out-of-phase vertical inputs to front and rear 
suspensions cannot be considered. 

The objectives of the present investigation were (1) to examine the relation between 
car spring and damping parameters and car responses when the car is subjected to simu- 
lated rail irregularities including sinusoidal and various transient and random inputs, and 
(2) to demonstrate a method of minimizing car accelerations by appropriate choices of 
spring and damping parameters for the vertical model. The minimization is more con- 
cerned with damping than with stiffness because damping properties are assumed to be 
more easily adjusted in an actual rail car than stiffness properties. 

The linear equations of motion were solved by digital computer programs for the 
sinusoidal inputs to both models and for one random input to the vertical model. To 
include nonlinearities in the lateral truck suspension system, a real-time simulation 
(RTS) digital computer program was developed and utilized. This program involved 
numerical integration of the equations of motion in scaled real time. Both linear and 
nonlinear acceleration responses to transient deterministic and random inputs were 
obtained for the lateral model by means of this program. 
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SYMBOLS 

Values are given in both SI and U.S. Customary Units. The measurements and 
calculations were made in U.S. Customary Units. 

Aq maximum acceleration response at steady state for qth degree of freedom 

a(t) generalized coordinate for vertical car bending 

B lateral distance between bolster springs (see fig. 13(b)) 

b(t) generalized coordinate for lateral car bending 

[cT] damping matrix 

Cq viscous damping constant for qth degree of freedom 

D dissipation energy 

d distance from end of car to center line of trucks (see figs. 1 and 13(a)) 

E Young r s modulus 

f forcing frequency, £- 

2 TT 

C0 Q 

fq frequency of qth degree of freedom, 

G modulus of rigidity 

g acceleration due to gravity 

g c structural damping coefficient of car in bending or torsion 

g rm s ro °t mean square of output spectral density in g units 

h-p vertical distance between car elastic axis and line of action of lateral trans- 

former spring (see fig. 13(a)) 
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vertical distance from traction-motor center of gravity to line of action of 
lateral equalizer spring (see fig. 13(b)) 

vertical distance from traction-motor center of gravity to line of action of 
lateral bolster spring (see fig. 13(b)) 

vertical distance between car elastic axis and line of action of lateral bolster 
spring (see fig. 13(b)) 

vertical distance of car section center of gravity from car elastic axis, posi- 
tive for center of gravity above elastic axis (see fig. 13(b)) 

flexural moment of inertia of car cross section 

mass moment of inertia of traction motor in roll about its center of gravity 
mass moment of inertia of car in roll about car elastic axis, per unit length 
torsional constant of car cross section 
stiffness matrix 

spring constant for qth degree of freedom 

length of car (see figs. 1 and 13(a)) 

mass matrix 

traction-motor mass 

car mass , mL 

transformer mass 

mass of car per unit length 

viscous damping coefficient for qth degree of freedom 


general displacement forcing function 


amplitude of forcing function 


Qo 

q(0 

<*2 - « 2 Q oJ 

R 

S(t) 

S^f) 

Sq(f) 

T 

t 

U 

v 

v(x,t) 

W(x) 

w(x,t) 

X 

Y(t) 

Yi(f) 

y A (t) 

y c (t) 


generalized coordinate 

amplitudes of sinusoidal motion in qth degree of freedom in terms of dimen- 
sionless quantities qj and 

lateral distance between equalizer springs (see fig. 13(b)) 
vertical displacement forcing function at rails (see fig. 1) 
input spectral density of vertical model (see eq. (3)) 
output spectral density of vertical model (see eqs. (7) and (8)) 
kinetic energy 
time 

potential energy 
car speed 

lateral bending deformation of car (see eq. (B2)) 

car bending-mode-shape value at any x 

vertical bending deformation of car (see eq. (A3)) 

length along car, measured from rear of car 

lateral displacement forcing function at rails (see fig. 13(b)) 

spectral density of lateral input at rails 

generalized coordinate of lateral translation of traction motor 
generalized coordinate of lateral rigid-body translation of car (see fig. 13(b)) 
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z A (t > 

z c (t) 

Z g (t) 

®m^m 
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Ay(t) 

6(t) 
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e A (t) 

0 c (x,t) 

X 

O) 

CO 

q 

Subscripts: 

A 

B 


generalized coordinate of lateral transformer translation 

generalized coordinate of vertical translation of traction motor 

generalized coordinate of vertical rigid -body translation of car (see fig. 1) 

generalized coordinate of vertical transformer translation 

eigenvalue properties of beam in bending (see eqs. (A8) and (A9)) 

center-to-center distance between axles of a truck suspension system (see 
eqs. (11) to (14)) 

displacement relation in lateral bolster springs for nonlinear spring behavior 
(see eqs. (15)) 

increment of vertical displacement at rails associated with cross -level (or 
rolling) displacement input (see fig. 13(b)) 

spectral density of cross -level input at rails 

generalized coordinate of car in rigid-body yaw 

generalized coordinate of traction- motor rolling degree of freedom 
rigid-body roll and torsion of car (see eq. (B3)) 
generalized coordinate for rigid-car roll 

wavelength for wave input 

angular forcing frequency, 2irf, rad/sec 

angular frequency of qth degree of freedom, 27rfq, rad /sec 


identifies properties associated with the traction-motor-equalizer system 
denotes vertical car bending frequency 
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associated with car 


identifies structural damping coefficient of car in lateral bending; also, 
lateral car acceleration (see table El) 

identifies structural damping coefficient of car in torsion 

identifies structural damping coefficient of car in vertical bending 

identifies forward trucks 

associated with transformer 

identifies energies and car bending stiffness of lateral model (see appendix B) 
identifies properties of traction-motor-lateral-equalizer system 
denotes lateral car bending frequency 

identifies properties of lateral bolster spring-damper system 
identifies properties of lateral transformer spring-mass system 
integer associated with mth car bending mode 
identifies amplitude of forcing function 

identifies bolster rigid-body damping coefficients and frequency 
identifies general degree of freedom 
identifies rear trucks 

identifies properties of vertical bolster spring-damper system 


associated with car torsion 



V 


identifies energies and car bending stiffness of vertical model 


a identifies roll stiffness and damping properties of traction motor 

Dots over quantities denote differentiation with respect to time. Bars over quanti- 
ties denote maximum amplitude in acceleration transfer functions. 

ANALYSIS OF VERTICAL MOTION 
Equations 

Mathematical model . - The mathematical model used in the analysis of the vertical 
motions of the railroad car is shown in figure 1. The model consists of a uniform beam 
representing the railroad car, a heavy mass Mg attached to the center of the car by a 
spring and damper (this mass represents an electrical transformer used in powering the 
railroad car), and two masses M^ attached by springs and dampers to the car and to a 
vibrating base representing the rails. The masses Ma represent the flexibly mounted 
sections of the truck frame and traction motors in the suspension system of the railroad 
car. For convenience in further discussions, each mass Ma with its attached springs 
and dampers is referred to as a truck suspension, the mass M^ is called the traction 
motor, the spring-damper combination between this mass and the car is called the 
bolster, and the spring-damper combination connecting the traction motor Ma to the 
rails is called the equalizer. A schematic view of the truck suspension system is shown 
in figure 2. The wheels, axles, and side frames (supporting the equalizer springs) are 
considered to be part of the rails, and the connections between the axles and traction 
motors are not taken into account in the model. In addition to the viscous dampers in the 
transformer and truck suspensions, structural damping is included in the car. The 
degrees of freedom considered are car rigid-body displacement, z c (t); car bending 
motion, w(x,t) = a(t)w(x); transformer displacement, Zg(t); and traction-motor displace- 
ment, z^(t). The model is symmetric about the middle of the car, and the same rail 
input is applied simultaneously to both truck suspensions; therefore, both traction motors 
have the same motion, and only one degree of freedom z a is required. With different 
inputs to front and rear trucks, separate degrees of freedom are required to specify the 
vertical motion of the front and rear traction motors, and rigid-body pitch would also 
have to be included. 

The equations of motion of the mathematical model are derived in appendix A. The 
matrix equations of motion are of the form: 

[M]{q> +[c](q> +[K]{q} = (q) + (q) <U 
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where the column matrix contains the displacements a(t), z c (t), Zg(t), and 

z A (t). In equation (1) [m] is the mass matrix, [c] is the damping matrix, [if] is 
the stiffness matrix, {q} is the displacement input, and {q} is the velocity input. 

Rail inputs .- The time history of the rail input used in this analysis was either 
sinusoidal or random. The sinusoidal input is 


Q(t) = S 0 sin cut 
Q(t) = S Q (x) COS cot | 


J 


( 2 ) 


where cu is the input frequency and S Q is the displacement amplitude. The random 
input is given by the spectral density variation identified as "middle estimate" in figure 8 
of reference 4 and is based on measured railroad roughness data. The spectral density 
for this displacement is 


Si(f ) = 


46.9 X 10-5 
f2.64 

7.27 X 10-5 
{2.64 


cm2/Hz 
in 2 /Hz ( 


(3) 


f being the input frequency in hertz. The root-mean-square (rms) amplitude of this 
spectral density is 1.113 mm (0.0438 in.) for a frequency range of 0.1 to 10 Hz. 


Acceleration transfer functions .- With sinusoidal input (eqs. (2)) and the assump- 
tion of sinusoidal motion in the form 


(q(t)J> = {q^ sin cut + cos cut 

equation (1) leads to the equations: 

[H - “ 2 M] - "0{5 2 } = {Qj 1 

> (4) 

“0(5.) + [M - » 2 [m]]{5 2 } =<Q 2 } J 

The elements of the matrices jjc], [m], [c], ^q-jJ>, > and {Q^} are 

given in appendix A. Equations (4) are solved for and {^ 2 } ^y using standard 

digital computer routines. From the solution of equations (4), the amplitude is 

determined by use of the relation 


(5) 


= /^lj + ^2j sin ( Wt + *j) 
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where Jq?. + q~~ is the response amplitude and xf/. = tan 1 is the phase angle. 

y t] 4i j qij 

In this paper, as in reference 4, no results are presented for the phase angle xpy 

On the basis of equation (5), a general expression for the steady-state acceleration 
response to sinusoidal input, or transfer function, may be written as 

Aqj = w 2 Q 0 + q| (6) 

The expressions for transfer function at different model locations are given in the fol- 
lowing table: 


Location 


Car 


Transformer 


Traction motors 


Ac(x) 

g s o 


Acceleration transfer function 


U)‘ 


g \ 


-|<s r-_ 

a.jW(x) + z c ^ + a2W(x) + 


A g oj 2 1-2 . -2 


gS 0 = g f gl + Z g2 


A to* -2 -2 

^ = tJ Z A1 + Z A2 



The two car locations for which acceleration transfer functions were studied were the end 
of the car (x = 0) and the middle of the car (x = L/2). 


Spectral density .- The spectral density of the acceleration response at a given 
model location is obtained from the transfer function and the spectral density of the 
input Sj(f) by the equation 



(7) 


The root-mean-square value of 


^rms - 



Sq(f) 
1 1/2 


in g units is obtained from 


(8) 


The integration limits f b and f e were taken as 0.1 and 10 Hz, respectively, for this 
study. 
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Optimization Procedure 

Sinusoidal input .- With the sinusoidal input, the objective of the optimization was to 
minimize the largest value of the acceleration transfer function occurring between 0 Hz 
and 20 Hz. The procedure used was the same as that used in reference 4. The first step 
in this procedure is to generate transfer function curves such as shown in figure 3 for a 
range of values of one of the system parameters. The next step is to plot the value of the 
transfer function associated with a given peak (determined from curves such as those of 
fig. 3) as a function of the varying system parameter. For example, figure 3(a) shows 
that peaks occur at about 1 Hz and 8.3 Hz in the transfer function at the end of the car, 
and figure 4(a) shows the variation of the transfer function associated with each of these 
peaks as the transformer frequency is varied. Similarly, parts (b), (c), and (d) of 
figures 3 and 4 show where the peaks occur for other positions in the model and how the 
peak values vary with transformer frequency. The frequency range associated with each 
peak is shown next to each curve in figure 4. 

The next step in the optimization is to identify a curve which shows the largest 
value of transfer function for each value of the varying parameter. Such a curve is 
called a maximum acceleration function and is indicated for each model location in 
figure 4 by the segments of the curves joining the circular symbols. A value of the 
varying parameter that minimized the maximum acceleration function is then chosen. 
Figure 4 shows that a transformer frequency of 6.5 Hz minimizes the maximum acceler- 
ation function at both car locations but not at either the transformer or traction- motor 
location. In the present study, a value was chosen that minimized the car location 
accelerations because they were considered more important from the standpoint of 
passenger comfort. 

Finally, with the optimum value of the first system parameter, a second system 
parameter is varied and the process described is used to determine an optimum value 
of the second parameter. In the present study, this procedure was repeated through four 
parameters and constituted an optimizing cycle. The cycle was repeated, as in refer- 
ence 4, until the value obtained as optimum for a given parameter was the same as the 
value obtained in the previous optimizing cycle. Fortunately, no more than two opti- 
mizing cycles were required in this study. Figures 4 to 7 show the curves used to deter- 
mine the optimum values of the varied parameters (with f^ = 4.12 Hz) and are discussed 
in the section "Results of Vertical Motion Analysis." 

Random input .- The objective of the optimization with the random input was to 
minimize the value of the root-mean-square (rms) acceleration in the frequency band 
from 0.1 to 10 Hz. The procedure was to obtain values of rms acceleration for each 
model location by using equations (4), (7), and (8) for a range of values of a system 


11 



■I 


parameter. A value of the system parameter was chosen to minimize accelerations at 
all locations, if possible, and the value chosen was then held constant while a second 
system parameter was varied. Only a single optimization cycle was required to obtain 
sufficiently accurate optimum values of the two parameters varied in this study. 

Results of Vertical Motion Analysis 

Two studies were carried out with the vertical model: (1) an optimization study 

based on minimum response of the car and (2) a parameter variation study of the effect 
on acceleration transfer function of two parameters that were not varied in the optimiza- 
tion study. The full scope of both of these studies is indicated in table I. Table 1(a) 
contains the values of the constant parameters and table 1(b) indicates the parameters 
that were varied, along with the value of each parameter used for each figure. 

Optimizati on for sinu soidal input.- The optimization procedure described for the 
sinusoidal input was used to determine optimum values of the frequency of the trans- 
former and the damping coefficients of the transformer, bolster, and equalizer. The 
curves used to determine optimum values are shown in figures 4 to 7, and the values of 
the remaining parameters are given in table I. The optimum transformer frequency was 
chosen from figure 4 to be 6 Hz. Figure 5 shows that the tranformer damping value 
of 0.25 minimizes the transfer function for both car locations and, in addition, results in 
near-minimum values of transfer function for the transformer and traction-motor loca- 
tions. Figure 6 shows that the transfer function is minimum at the end of the car for a 
bolster damping coefficient n s of 0.09 but at the middle of the car for n s = 0.14. The 
transfer function is sensitive to variations in bolster damping. In this case, n s was 
chosen as 0.11. This value is good for the transformer location but not for the traction 
motor, for which the minimum transfer function is obtained for n g = 0.5. Figure 7 
shows that optimum equalizer damping is 0.03 for the end of the car but is 0.31 for the 
middle of the car. The compromise value of n^ = 0. 1 was chosen arbitrarily, between 
the two minimizing values. 

Figures 4 to 7 show that it was not possible to choose a value of the varying param- 
eter that minimized the transfer function at all model locations simultaneously. This 
result was also observed for the three-degree-of-freedom model of vertical car motion 
studied in reference 4. Because of this result, care must be taken in a real rail car 
situation that an optimum value of a system parameter chosen to minimize acceleration 
at, for example, a car location does not result in excessive acceleration at another loca- 
tion, such as the traction motors. The factors to be considered when the accelerations 
at both car locations cannot be minimized together, such as those shown in figures 6 
and 7, are discussed in reference 4. 
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The results shown in figures 4 to 7 were obtained by using an equalizer spring 
value that resulted in an equalizer frequency of 4.12 Hz. The study represented by 
these figures was repeated with an equalizer frequency of 5.53 Hz (a stiff er equalizer 
spring). The values of the other nonvarying parameters were the same as those shown 
in table 1(a) for figures 4 to 7. The optimizing curves with the equalizer frequency 
fA = 5.53 Hz were qualitatively very similar to the curves shown in figures 4 to 7; the 
optimum values of the varied parameters are compared for the two equalizer frequencies 
in table II. This table shows that the optimum values of only transformer and bolster 
damping were different for the two equalizer frequency values. Moreover, the transfer 
function with = 5.53 Hz was insensitive enough to variations of transformer 
damping ng to insure that a value of rig = 0.25 would be an acceptable value. 

Figure 6 shows, however, that transfer function is very sensitive to variations of bolster 
damping ng; therefore, the difference between n s = 0.057 and n s = 0.11 is signifi- 
cant. These results indicate (as might have been expected) that the value of the equalizer 
spring is important in determining one of the optimum damping values, bolster damping. 
In this particular case, the equalizer spring value is not important for the other param- 
eters fg, ng, and n^. 

Optimization f or random input .- By using the optimization procedure previously 
described for the random input, optimum values of bolster and equalizer damping coef- 
ficients were determined. The results are presented in figures 8 and 9. The values of 
all other parameters were the same as those for the sinusoidal input study of figures 4 
to 7, as shown in table I. The values of transformer frequency and damping were not 
varied; however, the optimum values from the sinusoidal input study were used. 

Figures 8 and 9 both show that the rms acceleration, in general, was insensitive to varia- 
tions of bolster and equalizer damping. Figure 8 shows that a bolster damping coeffi- 
cient of about 0.3 is a good value for minimizing rms accelerations. In spite of the 
insensitivity of rms acceleration to variations of bolster damping, the value of n s = 0.11 
as determined from the sinusoidal study has rms accelerations at the car location about 
30 percent greater than those for n s = 0.3. Therefore, these rms studies indicate that 
a higher bolster damping value is required to minimize rms accelerations than is 
required to minimize peak acceleration transfer functions. This result is in agreement 
with the results for the three-degree-of -freedom model of reference 4, where, with no 
mass between the rails and the car, it was shown that the optimum suspension damping 
values were larger for random inputs than for sinusoidal inputs. Figure 9 shows, 
however, that rms acceleration is so insensitive to variations of equalizer damping (for 
the particular values used for the other system parameters) that any value between 
n A = 0 and nA = 0.5 may be chosen without penalty to rms acceleration. 
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The rms values shown in figures 8 and 9 were obtained from spectral densities of 
the type shown in figure 10. These particular spectral densities apply to the optimum 
values of ns and n^ = 0.30 from figure 9 and show a concentration of energy at the 
low frequency range (around 1 Hz) for the two car locations and the transformer. For 
the traction motors, the energy concentration occurred in the higher frequency range, 
near 5 Hz. 

Parameter variat ion.- Car bending frequency fg and car weight M c g were held 
fixed at 8 Hz and 574 kN (129 000 lb), respectively, for the optimization studies, and the 
structural damping coefficient g cv was taken equal to zero. To determine how these 
properties might interact with other parameters in the model, including those chosen 
for optimization, transfer functions were obtained for wide variations of car bending 
frequency and car weight and for a nonzero structural damping coefficient. Maximum 
values of the peak acceleration curves (maximum acceleration functions as described in 
section "Optimization Procedure") and their associated frequencies are plotted for the 
two car locations as functions of bending frequency in figure 11 for an arbitrary non- 
optimum configuration. Variations of acceleration transfer functions with car weight are 
shown in figure 12 and were obtained by varying both car weight and bending frequency, 
car bending stiffness being assumed to be fixed at a value corresponding to fg = 8.0 Hz 
and M c g = 574 kN (12S000 lb). The values of all parameters used to generate fig- 
ures 11 and 12 are shown in table 1(b). 

Figure 11 shows that the acceleration transfer function with no car damping has a 
maximum value for both car locations for a car bending frequency value of 3 Hz. The 
value of transfer function at its maximum is about 8.5 and 3.9 times the value at 
frequencies above 6 Hz for the end of the car and the middle of the car, respectively. 

This sensitivity suggests the need for knowing the bending frequency of a given rail car 
and its possible variations in service so that large values of acceleration can be avoided 
by proper design of the system. Figure 12 shows that car acceleration transfer functions 
were relatively insensitive to wide variations in car weight. 

Comparison in figure 11 of the solid curves for zero structural damping with the 
dashed curve for structural damping ^g cy = 0.05^ shows that* a small amount of structural 
damping causes large decreases of acceleration transfer function within the frequency 
range of resonance (<6 Hz). However, if, by proper design, this range can be avoided, 
structural damping becomes of minor importance, as indicated by the virtual coalescence 
of the curves of figure 11 away from resonance. 
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ANALYSIS OF LATERAL MOTION 


Mathematical Model 


Views of the mathematical lateral model and its coordinate system are shown in 
figure 13. Viscous damping representation is omitted for figure clarity. The motion of 
the model is governed by equation (1), and the matrix elements for this equation are 
derived in appendix B in the same manner as that used in appendix A for the vertical 
model except for differences in certain strain-energy terms due to the nonlinear lateral 
bolster springs. Car motion is represented by the torsional and lateral bending of a 
uniform beam. As in the case of the vertical model, linear viscous damping is assumed 
for the lateral truck and transformer suspensions and structural damping for the car 
itself. Symmetry about the center of the car is assumed in lateral as well as longitu- 
dinal directions. 


Degr ees of freedom .- The lateral model consists of ten degrees of freedom identi- 
fied as follows: 

Car lateral bending mode, b(t) 

Car rigid -body lateral translation, y c (t) 

Car rigid -body roll, 0 C q (t) 

Car first torsion mode, 

Car rigid -body yaw, r]( t) 

Transformer lateral translation, yg(t) 

Lateral translation of rear traction motor, y^ r (t) 

Lateral translation of forward traction motor, y^(t) 

Roll of rear traction motor, «Ar« 

Roll of forward traction motor, 0^(t) 


The car is free to twist about a longitudinal elastic axis (designated in fig. 13(b)), and 


the first torsion mode shape of the car is given by cos^j. The total bolster spring 
constant k s and damping constant c s for each truck are divided equally between the 
springs and dampers located symmetrically beneath each side of the car, as shown in 
figure 13(b). Similarly, the total equalizer spring constant k^ and damping con- 
stant c^ in each truck are each divided equally between the springs and dampers 
connecting the traction motors to the rails. Both rigid-body roll and torsion M kA 
and 0 c i(t)) are resisted by the bolster springs and dampers. The spring constant 
and damping constant over each rail represent the combined stiffness and damping 

of separate equalizers actually located near each wheel; that is, there are four equalizers 
per truck. Traction-motor roll ( 0 Ar (t) and ft)) * s assumed to be resisted by this 
spring-damper system. Additional stiffness and damping in the equalizer -traction- motor 
system are provided by a rolling spring constant k ff and a damping constant c a at 
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the traction -motor center of gravity, as shown in figure 13(b). The introduction of this 
additional stiffness and damping is intended to simulate the complicated linkage and gear 
mechanisms between the traction-motor and wheel-axle system in an actual rail car 
truck, as illustrated in figure 2. 

Lateral stiffness and damping in the transformer and truck suspensions are given 
by spring and damping constants that are assumed to approximate lateral shear stiffness 
and damping in the vertical transformer suspension, the bolsters, and the equalizers. 

The lateral transformer spring and damping constants are designated by k^ g and cj g , 
the lateral spring and damping constants of the bolster by k^ c and c i c , and the lateral 
equalizer spring and damping constants are denoted by k^ and c i . The transformer 
properties are assumed to be linear, and the bolster and equalizer properties are also 
assumed to be linear between the small solid semicircles in the truck suspension view 
(fig. 13(b)). These semicircles represent compressible stops or bumpers, each having 
the nonlinear spring function shown in figure 14. This function is based on measured 
force-deflection data for an actual truck suspension. 


Rail inputs .- The lateral model was excited by lateral and rolling, or rocking, inputs 
from the rails as indicated by the horizontal and vertical arrows at the bottom of the 
truck suspension view in figure 13(b) . The lateral input is designated Yf(t) and Y r (t) 
at the forward and rear trucks, respectively. The roll input is identified as cross level, 
as in reference 1, and was generated by the differences (t) and 5 r (t) between 
unequal vertical inputs applied to each rail as shown in figure 13(b). This difference was 
always applied as the displacement of the right-hand rail (with the car viewed from the 
front). In general, these inputs may be written in vector form as 
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Three classes of input were applied to the lateral model: sinusoidal, transient, 
and random as indicated in figure 15. The sinusoidal inputs (fig. 15(a)) may be written in 
the general form: 

Qf (t) = Q 0 f sin cot ^ 

Q r (t) = Q or sin ( wt + «) I (1Q) 

Qf(t) = Q 0 f w cos 

Q r (t) = Qor 00 cos + K b 

where co is the input frequency and Q^, Q or = Y 0 f, Y or or 6 0 f, 6 0r depending 
on whether the model is excited by lateral or cross -level inputs, respectively. The 
quantity k represents the phase angle between forward and rear inputs, and the model 
was excited with inputs either in phase or out of phase as illustrated in sketch (a). The 
out-of-phase option allowed for the influence of car torsion and traction-motor rolling 
effects to be examined. 
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The three nonsinusoidal deterministic inputs in figure 15 are transients repre- 
senting different types of disturbances of finite length X along the rails. The first of 
these inputs is a single square pulse (fig. 15(b)) with a constant displacement applied to 
the front trucks Q Q f and rear trucks Q or in a sequence of time intervals which intro- 
duces car speed v as an additional variable through the simple relations 


A + (AZ)x 

tl - t 0 = t 3 - 1 2 = v 

(11a) 

L - 2d - (AZL - X 


t 2 -ti= y 

(lib) 


where t^ - tQ is the time interval required for the forward truck to roll over the dis- 
turbance, t 2 - t^ is the time between when the front truck leaves the disturbance and 
the rear truck encounters it, and t 3 - t2 is the time it takes for the rear truck to 
roll over the distrubance. If L— - 2d - (AZ) t is less than X, corresponding to a long 
disturbance, t2 - t^ becomes negative and indicates that the rear truck encounters the 
disturbance before the front truck leaves it. The distance (AZ)^ is the center-to-center 
distance between the axles of a truck. The velocity inputs Qf(t) and Q r (t) are 
assumed to be zero. Lateral or cross-level inputs were applied in the directions indi- 
cated in the foregoing in-phase sketch, and the only phase lag of the rear truck behind 
the forward truck was due to car speed, in accordance with equations (11). 

The second transient input is the square pulse train (fig. 15(c)) which is intended to 
reproduce the effect of rail joints as illustrated in sketch (b): 


f— 

5 — 


* 




L - 2d -(AZ) t 


A — 


V 2 


h 




.Qli 


OO 


Sketch (b) 

The length Xr is the distance between joints, and X is the gap between the ends of the 
rail joints. In the present paper, X was chosen as 2.54 cm (1 in.) and X R as 
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11.9 m (468 in. or 39 ft). This input is treated as a succession of square pulses applied 
to one rail (the right-hand one for the front view of the car, consistent with the foregoing 
in-phase sketch). The displacement inputs to each truck are constant and the velocity 
inputs are assumed to be zero, as for the square pulse. In addition to equations (11), 
two other time-distance relations are involved, namely 

X + (AZ). 

t 5 - t 4 = (12a) 


and 



(12b) 


where t 5 - t 4 is the time required for the forward truck to roll over the disturbance 
due to the gap in the other (left hand) rail, and t 5 - t^ is the time required for a truck 
to traverse the distance between gaps which occur on opposite rails. 


The third transient input is a single-wave pulse (fig. 15(d)) which is intended to 
simulate the effects of a lateral deviation of the track from a straight course or the bank 
on a curve, depending on whether the inputs are lateral or cross level. This input is 
governed by trigonometric relations of a form similar to those used in reference 2. 

For the forward truck, the displacement input is 




cos 


2jrv(t - t 0 )' 
X + (AZ)j- 


(13a) 


and for the rear truck 


Q r (t) = 


Qpr 

2 



cos 


277v(t - t 2 ) 
X + (AZ) t 


(13b) 


The velocity input for this case is 


Qf(t) = 


Qof77V 

X + (AZ)t 


sin 


2frv(t - t Q ) 
X + (AZ) t 


(14a) 


Q r (t) = 


Qpr 7 ^ 

A. + (AZ)j- 


sin 


2ttv (t - t 2 ) 
X + (AZ) t 


(14b) 


Equations (11) apply for this input as well as for the square pulse. In the application 
of equations (13) and (14), X was varied over a wide enough range to insure that 
L - 2d - (AZ)j. - X was both positive and negative. 
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Random inputs to the lateral model included the functions illustrated in figures 15(e) 
to 15(g). Spectral densities identified as 1 and 2 in figure 15(e) correspond to rms ampli- 
tudes of 5.08 cm (2 in.) and 2.54 cm (1 in.), respectively. Spectral density 3 is based on 
an experimental measurement of random vertical displacements along a 3.2-km (2-mile) 
section of test track in the eastern United States, and the rms amplitude was 1.524 mm 
(0.06 in.) for a vehicle speed of 49.2 m/sec (110 mph). This input was arbitrarily applied 
to the lateral model with a 2.54-cm (1 in.) rms lateral amplitude. Figure 15(f) illus- 
trates the square-pulse train input with random amplitudes. Figure 15(g) shows a time 
history of a combination of this random-pulse train input with spectral density 4 in fig- 
ure 15(e) for a 2.54-cm (1 in.) rms amplitude. This input was also applied arbitrarily 
as lateral input to the lateral model. Spectral density and rms inputs are denoted 
by Yj(f) and Y rms for the lateral input and 6j(f) and 6 rms for the cross-level 
input. All these random inputs were applied simultaneously to front and rear trucks. 
Velocity random inputs were neglected. 

Outputs . - Responses of the lateral model were obtained in the forms of acceleration 
transfer functions for the sinusoidal inputs, absolute values of maximum accelerations 
in g units for the transient inputs, and acceleration spectral densities and rms values 
for the random inputs. Lateral acceleration responses were calculated in the car at 
various stations distributed along the longitudinal center line of the car and also for the 
transformer and traction motors. Vertical acceleration responses were determined in 
the car just over the bolsters and in the traction motors over the equalizers. Accelera- 
tion transfer functions for the sinusoidal input were obtained from equations (4) to (6) as 
for the vertical model, and particular forms of equation (6) are listed in table HI for the 
lateral model. 

Solutions of equations of motion.- Equations (4) and the transfer functions listed in 
table HI were programed on the Control Data 6600 computer system at the Langley 
Research Center, as were the corresponding equations for the vertical model. This 
program is restricted to linear spring elements. The solutions of equation (1) for the 
transient and random inputs in figure 15 were programed, along with equations (11) to 
(14), for another Langley computer facility, discussed in the following section, that was 
also able to handle the nonlinear spring characteristics shown in figure 14. 

Scaled Real-Time Simulation 

The presence of the nonlinear springs and the wide selection of inputs previously 
described prompted application of another method of solving equation (1) for the lateral 
model. This method involved the numerical integration of this equation in time by use of 
the real-time simulation (RTS) capabilities of the Langley central digital data recording 
facility. In this approach, the ten simultaneous equations represented by equation (1) 
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were integrated in the time domain on the Control Data 6600 computer system, and solu- 
tions were fed through digital-to-analog converters (DAC) to give time histories of 
displacements, velocities, and accelerations in various parts of the model. To insure 
sufficiently accurate integrations, the independent variable time had to be slowed, or 
scaled, during the computation. A time interval of 1.5625 msec was used in a second- 
order Runge-Kutta integration scheme. Results from this program for sinusoidal 
inputs were checked with those obtained from the digital transfer -function program. The 
sinusoidal inputs of equations (10) could be applied in the RTS program for an arbitrary 
phase angle /c, whereas the digital solution was programed only for solutions with k 
equal to integral multiples of tt corresponding to in-phase and out-of-phase inputs to 
the front and rear trucks, respectively. 

An important feature of the RTS approach is the on-line random input capability 
represented in the flow chart of figure 16. In this arrangement, band-limited white 
noise, which was generated by an analog computer, could be introduced directly through 
a cut-off filter, or through a desired shaping filter, into the RTS computer program for 
solving equation (1). Outputs of this program and the random input were both fed on-line 
into a time series analysis digital program (TSA), in which statistical properties of input 
and output were calculated. Spectral density and rms value were the two properties of 
interest in this investigation. By appropriate choices of shaping filters, reasonable 
approximations to the spectral-density shapes for inputs 3 and 4 in figure 15(e) could be 
generated for any desired input amplitude. 

Real-time simulation system .- The real-time simulation (RTS) system used to 
solve the ten equations of motion provides console control combined with digital accuracy 
and convenience. (See ref. 6.) The ten second-order derivatives in equation (1) were 
formed and numerically integrated twice during each machine iteration to yield the 
velocities and displacements. These operations are indicated on the flow chart of 
figure 17. During most of this study, 32 iterations per second (20-to-l slow) was the 
normal iteration rate, whereas 64 iterations per second (10-to-l slow) were used when- 
ever possible. 

Ran dom studies .- An Electronic Associated, Incorporated (EAI) 231-R analog 
computer was programed with appropriate filters and a mean and rms control circuit to 
produce desired random inputs to the mathematical model from an analog Gaussian 
noise generator. This equipment is represented by the random-input (or noise gener- 
ator) box in the flow chart of figure 16. Outputs from the mathematical model were 
directed to an on-line TSA program located at a separate control point of the RTS system. 
This program computes mean values, standard deviations, chi-square tests for normality, 
histograms, autocorrelation functions, and power spectral density functions (PSD) for its 
input channels. These statistical properties are discussed in detail in reference 7. A 
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CDC-250 cathode ray tube (CRT) was used to display linear or log -log plots of the PSD 
values generated by the TSA program. Control stations such as those used for solving 
the equations of motion, operation of the TSA program, and operation of the CRT are 
shown in figure 18. Also shown are a CRT console and the time-history recorders, each 
capable of displaying 8 different displacements, velocities, or accelerations. 

Visual car display.- An oscilloscope display showing a front cutaway view of the 
car at the forward truck was driven by a program on the GPS 10 000 iterative analog 
computer. As indicated in figure 16, this analog program transformed DAC outputs from 
the mathematical-model solution into signals for the scope display. This display was of 
small analytical value because of the 20-to-l or 10-to-l time scales but was useful for 
demonstration purposes. 


Results of Lateral Motion Analysis 

The presence of eight spring elements and eight dampers in the lateral model 
requires the use of a computerized algorithm for minimizing vibration levels in the 
system. The trial -and- error cycling approach utilized in a limited optimization of the 
vertical model, with eight system elements at most, was considered to be impractical 
for the lateral model. Instead, the analytical tools previously described were exploited 
to obtain trends of acceleration responses to inputs involving realistic ingredients for 
lateral and cross-level rail disturbances. Parameter variations included realistic rail 
car values, and the widest variations were made for the sinusoidal inputs utilizing the 
transfer-function digital program. The RTS capabilities were employed to explore 

(1) relations between car speed and damping coefficients in the truck suspension and 

(2) nonlinear spring effects for the transient and random inputs. Some effects of input 
wavelength were also considered. The full scope of this study is given in tables IV 
and V. The fixed properties are listed in table TV, and the values of the variable prop- 
erties and inputs in each of the figures are presented in table V. Inputs in table V are 
identified according to the appropriate parts of figure 15. 

Results for sinusoidal input.- Typical acceleration transfer functions for various 
stations in the lateral model are shown in figure 19. With all damping coefficients 
reduced to 0.001 and the car structural damping coefficients set equal to zero, peak 
acceleration responses occurred near the uncoupled frequencies of the system. The 
existence of only a few peaks in figure 19 is attributed to the moderate to high damping 
levels indicated by the values listed in table V for this figure. A dominant peak 
at 0.44 Hz occurred frequently in the car and transformer responses as shown in 
figures 19(a) and 19(b). The response was sensitive to the magnitude of the lateral 
bolster damping coefficient nj c . Other parts of the model experienced dominant peak 
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accelerations at higher frequencies for the cases considered, as is evident, for example, 
in figures 19(c) to 19(e). 

All ten degrees of freedom were retained in nearly all the studies; however, the 
effects of deleting traction-motor rolling degrees of freedom (0^ r (t) and 0^(t)) and of 
changing traction-motor stiffnesses (k^j. and were given limited attention. The 

effects of these parameter variations are compared in the transfer -function plots of 
figures 20 and 21. Figures 20(a) and 20(b) show only a 7-percent reduction in maximum 
lateral car accelerations to lateral in-phase inputs because of the complete elimination 
of these rolling degrees of freedom and about 4-percent reduction with traction-motor 
rolling stiffnesses increased from k^ = k a f = 26.8 to 101.5 MN-m/rad (corresponding 
to 2.375 x 10® to 9.00 X 10® in-lb/rad and to a frequency increase from fQ, = 18 Hz to 
35 Hz). Somewhat larger effects on the vertical bolster responses are evident in fig- 
ure 20(e). Figure 21 shows still larger effects of traction-motor rolling stiffness on 
maximum acceleration responses to cross-level out-of -phase inputs. Lateral end-of- 
car response (fig. 21(a)) and the vertical response over the bolster (fig. 21(c)) are each 
reduced by about two-thirds with f a increased from 18 to 35 Hz; however, these 
maximum responses occurred close to 18 Hz, and at lower frequencies within the range 
of interest of this study, these effects are seen to be much smaller. In general, changes 
in car accelerations due to these variations in traction- motor stiffness and rolling 
degrees of freedom in figures 20 and 21 are considered to be small by comparison with 
the effects of other parameters in the model. 

Maximum values of peak acceleration curves are shown as functions of frequency 
(that is, spring) and damping parameters in figures 22 to 29 for sinusoidal inputs. 
Frequencies associated with these peak accelerations are identified at points along the 
curves. Figure 22 reveals a resonant region of significantly increased lateral car 
accelerations due to lateral car bending frequencies (eq. (B16a)) falling within range of 
other frequencies of the model. These trends are similar to those observed in figure 11 
because of vertical bending frequencies. The great reduction in lateral car accelerations 
due to car structural damping in this resonant region is similar to that observed in 
figure 11 for the vertical model. Likewise, with lateral car bending frequencies well 
removed from the other frequencies of the system, lateral car accelerations were 
greatly reduced and were insensitive to variations in structural damping, as were verti- 
cal car accelerations in figure 11. For a lateral bending frequency of 8 Hz, the maxi- 
mum transfer function was about 0.25 g per inch at both the end and middle of the car and 
occurred at a frequency of 0.44 Hz. The results in figure 22 were obtained with lateral 
in-phase inputs from the rails. 

Maximum lateral car acceleration responses to out-of -phase lateral and cross- 
level inputs are plotted as functions of car torsional frequency in figure 23. Lateral 
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responses in figure 23(a) are seen to be constant over a wide range of torsional fre- 
quencies and significantly reduced with fx less than about 5 Hz. The effect of cross- 
level inputs in figure 23(b) is somewhat reversed, maximum acceleration responses 
increasing as fx is moved beyond the range of other system frequencies. It may also 
be observed that the lateral car responses to lateral inputs were higher than the lateral 
car responses to cross-level inputs, and this result is a general one of this study not 
only for the sinusoidal inputs but for the other inputs as well. Higher car responses to 
lateral than to cross-level inputs were also obtained in reference 1. The large varia- 
tions in response shown in figures 22 and 23 indicate the need for accounting for car 
lateral bending and torsional stiffnesses in analyzing car vibrations. Comparison of the 
trends in figure 23, particularly figure 23(b), with those in figure 22 suggests that tor- 
sional stiffness plays a smaller role than bending stiffnesses in car vibrations. 

Variations in acceleration transfer functions with car weight are shown in figure 24 
for out-of-phase inputs. These variations were obtained by varying car weight, rolling 
moment of inertia, lateral bending frequency, and torsional frequency with car bending 
and torsional stiffness assumed to be fixed at values corresponding to f^g = 8.0 Hz, 
f T = 15.0 Hz, Meg = 574 kN (129000 lb), and Ig^ = 3090 kg-m (695 lb-sec2). Although 
the trends differ from those shown for the vertical model in figure 12, the effects of a 
wide variation in car weight on lateral car transfer functions are small, as these effects 
were for the vertical car transfer functions. 

The flexibly mounted transformer beneath the middle of the car acted as a vibration 
absorber for reducing lateral car accelerations in an analogous manner to its role in 
reducing vertical car accelerations in the vertical model. The peak responses in 
figures 25 and 26 show reductions that were obtained for a wide range of transformer 
frequencies and damping coefficients with all other parameters held constant at the 
values listed in tables IV and V. 

A vibration-absorber role can also be assumed for the lateral truck suspension. 
However, the present trend study was limited to the effects of damping with fixed spring 
parameters. Thus, with fixed lateral bolster and equalizer frequencies, lateral bolster 
and equalizer damping coefficients were varied and reductions in peak acceleration 
curves obtained as in figures 27 and 28. The results in figure 27 show relative minimum 
lateral car accelerations at different car locations for different bolster damping coeffi- 
cients due to lateral in- and out-of-phase inputs. In this case, structural damping has 
the effect of widening these differences in the bolster damping values for minimum 
acceleration responses. Similar results are shown in figure 28 for the equalizer damping 
variation, except for the effect of structural damping which is not included. Figure 29 
shows the relative minimum of peak lateral car accelerations at different values of 
bolster roll damping coefficient for in-phase and out-of-phase cross -level inputs. With 
car torsional structural damping included, out-of-phase responses were sharply reduced 


24 



with a minimum lateral response indicated at n 0 T = 0. This reduction in response due 
to structural damping is consistent with that shown in figure 23. These results indicate 
the joint influence of both car and suspension damping on lateral car acceleration 
responses. 

Results for transient a nd random inputs .- In utilizing RTS capabilities to examine 
linear car responses to various transient deterministic and random inputs, the main 
variables of interest were car speed v, input disturbance length X, and lateral bolster 
damping coefficient nj c . Some effects of bolster roll damping coefficient n Q <p were 
also explored. Most of the results are presented in plots of peak lateral car acceleration 
responses as functions of car speed, as shown in figures 30 to 35. Some variations of 
peak responses with car length are shown in figure 36 for the square-pulse train input of 
figure 15(c). 

These results show wide differences in peak lateral car responses due to changes 
in input and changes in car speed for individual inputs. Lateral inputs generally produced 
higher lateral car responses than did cross-level inputs, as was noted in figure 23 for 
sinusoidal inputs. Mid-car responses were lower than end-of-car responses for tran- 
sient deterministic inputs (figs. 30 to 33) and for the combined random input (fig. 35). In 
some cases, such as those illustrated in figure 36, minimum car responses along the 
length of the car occurred at about 25 percent and 75 percent of the car length. More- 
over, the spread between front and rear responses of the car was considerably larger for 
the wave-pulse inputs than for the other inputs. Figures 32 and 33 show responses at the 
rear of the car that are significantly higher than those at the front of the car. Sample 
time histories in figures 30 and 32 further illustrate differences in the character of the 
car responses for square-pulse and wave-pulse inputs for two widely different distur- 
bance lengths. Initial accelerations due to a pulse input for X = 0.305 m (12 in.) in fig- 
ure 30 were more erratic than those due to a wave input for X = 76.2 m (3000 in.) in 
figure 32. Peak accelerations were about 0.047g for both time histories. 

Other indications of how dependent these results were on the nature of the input are 
afforded in comparisons of rms accelerations in table VI for two different random inputs, 
and in comparisons of acceleration responses for deterministic and random pulse 
train inputs in figure 34. The table compares rms lateral responses for the test track 
input, corresponding to input spectral density 3 in figure 15(e), and the combined input of 
figure 15(g) for a 2.54 cm (1 in.) rms amplitude. The results are taken from figures 35 
and 37 for the same car speed 49.3 m/sec (110 mph) and lateral bolster damping coeffi- 
cient ^nj c = 0.11). Responses to the test track input are seen to be 22.5 to 39 percent 
higher than responses to the combined input, and mid-car responses are higher than 
responses at the truck locations in contrast to the combined -input responses. In 
figure 34, n^ c was 0.057, and the lateral-input amplitudes were 2.54 cm (1 in.) for the 
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deterministic pulse train input and 2.54 cm (1 in.) rms for the random pulse train input. 
As the figure shows, lateral responses due to random input were larger than those due to 
deterministic transient input, and differences in lateral responses between the trucks and 
the middle of the car were greater for the random than for the deterministic transient 
inputs. The lateral responses to random pulse train inputs (X = 2.54 cm (1.0 in.)) 
followed about the same trends and were about the same magnitudes as the lateral 
responses to the deterministic pulse inputs (X = 0.305 m (12 in.)) in figure 30. However, 
rms accelerations at rear and forward truck locations tend to be more nearly equal in 
figure 34 than do the corresponding peak accelerations in figure 30; thus, the effects of 
input variation on car response are also apparent in this comparison. These results 
further support the recommendations in references 1 and 4 which point out the need for 
obtaining measured input spectra for studies of this kind. 

Effects of wavelength X were explored for the wave-pulse inputs, and figure 33 
shows the variation of peak lateral car accelerations with car speed for two different 
values of X, both less than that of figure 32. Comparison of the results in these two 
figures indicates that lowering X from 76.2 m (3000 in.) to 19.05 m (750 in.) substan- 
tially reduced the speed range bounding the wide spread in end-of-car acceleration 
responses to lateral inputs. The much smaller lateral response to cross-level inputs 
behaved in a manner somewhat similar to that of the lateral-input responses; however, 
the spread in end responses still extended over about the same speed range regardless 
of X value, and cross-level responses for X = 38.1 m (1500 in.) are shown in 
figure 33 to be nearly indistinguishable from those for X = 19.05 m (750 in.). 

In addition to the effects of input variations and disturbance length, variation of 
damping in the lateral bolster ( n £c) a ^ s0 caused large fluctuations in lateral car accelera- 
tion response. For example, localized increases in response were observed at particular 
car speeds for the deterministic pulse train and combined random inputs in figures 31 
and 35; however, the effects of this damping are reversed. That is, increased damping 
for the deterministic input tended to increase responses near 53.7 m/sec (120 mph, 
fig. 31), whereas increased damping for the random input not only gave reduced responses 
but also caused them to be largely invariant over the speed, range (fig. 35). It should be 
noted that car speed variations for the combined random input of figure 35 are due solely 
to the pulse train part of the input. By contrast with these effects, the wave-pulse input 
caused steep increases in lateral car peak responses to maximum values near 44.7 m/sec 
(100 mph) and slight decreases at higher speeds (fig. 32). Within the 35.8 to 44.7 m/sec 
(80 to 100 mph) speed range, increased lateral bolster damping coefficients were more 
effective in reducing lateral car acceleration responses than at other speeds. These 
damping effects clearly indicate that no fixed damping value would reduce lateral car 
accelerations at all speeds or at all locations in the car for any of the inputs considered. 
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Effects of bolster roll damping on lateral car acceleration responses were much 
smaller for cross-level inputs than the effects of lateral bolster damping were for lateral 
inputs. An indication of these effects may be seen in figure 30 where responses for 
n oT = 0-30 are higher than those for n 0 >p = 0.20 at nearly all car responses from 
17.88 to 71.52 m/sec (40 to 160 mph). Other results (not included) show still higher 
responses for n oT = 0.40. As observed for the lateral bolster damping effect, the effect 
of bolster roll damping was also very much a function of car speed. 

Nonlinear spring behavior .- All the results presented thus far are based on linear 
spring characteristics in the model. In addition to its application to the simulation of 
various rail irregularity inputs, the RTS method was also used to explore the effects of 
spring nonlinearity in the lateral bolster suspensions. The nonlinear spring characteris- 
tics shown in figure 14 were assumed for this purpose, and the deflections in the 
abscissa of the figure are related to the degrees of freedom in the model through the 
equations* 


Ay r (t) = -b(t) W(d) - y c (t) + (— - d) r]( t) + y Ar (t) - h 2r 0 Ar (t) 


- hi 


ird 


e c0 (t) + % ! (t) cos — 


(15a) 


Ayf(t) = -b(t) w(l - d) - y c (t) - (y - d ) v(t) + y Af (t) - h 2f e Af (t) 


- hi 


0 cO (t) + 0 cl (t) cos — 


(15b) 


Through these equations, nonlinear behavior in the bumpers directly below the car is felt 
in the car and the traction-motor suspension. 

With the nonlinear spring characteristics allowed for in equation (1), acceleration 
responses to sinusoidal inputs failed to reach a steady-state condition but varied in 
amplitude for a finite number of cycles of oscillation in repeated sequences in time. This 
response behavior ceased after termination of the sinusoidal forcing function. With band- 
limited lateral white-noise amplitude of 2.54 cm (1 in.) rms input (no. 2 in fig. 15(e)) 


Additional relations are included in reference 5 for an assumed nonlinearity in the 
lateral equalizer springs. However, these relations (eqs. (6) of ref. 5) have been found 
to be incorrect because of the omission of lateral inputs from the right sides of the equa- 
tions. Consequently, the "nonlinear" curves in figure 11 of reference 5 are incorrect and 
have been omitted from figure 37 of the present paper. In addition, the terms -l^r^ArW 
and -h2f0 A f(t) in equations (15) were erroneously omitted from equations (5) of refer- 
ence 5 but were included in calculations for the nonlinear spectral density in figure 7 of 
this reference. 
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to equation (1), output spectral densities and rms values, such as those shown in figure 38, 
were obtained and are seen to be considerably higher for the nonlinear case than for the 
linear case. Figure 38 also shows large reductions in responses due to increasing 
lateral bolster damping from nj c = 0.057 to 0.11. 

Similar effects of the nonlinear springs were observed in lateral acceleration 
responses along the length of the car due to wave pulse inputs as shown in figure 39. 

Here the input wavelength X was 9.52 m (375 in.), and the car speed 35.8 m/sec 
(80 mph). The solid curve identified as "linear” applies to lateral acceleration responses 
for which k i c was linear and is included for comparison with the other curves for 
which kj c was nonlinear. As may be seen, increasing n^ c from 0.057 to 0.30 
resulted in drastically reduced responses, as in figure 38 for the random input. With 
n i c = 0.30, the responses were reduced to the levels of the linear responses over the 
length of the car. 

It thus appears that higher car accelerations can be obtained with nonlinear springs 
than with linear springs in the model and that this result is not likely to change when the 
input changes. It is also evident that a moderate amount of damping can be very effective 
in reducing nonlinear acceleration responses. 

CONCLUDING REMARKS 

Simplified vertical and lateral dynamic models of a high-speed rail passenger car 
have been developed and studied to obtain dynamic response information which is relevant 
to the problem of passenger comfort. Vertical responses of the car to vertical inputs 
from the rails have been obtained in a four-degree-of-freedom model, and lateral 
responses to lateral and rolling (cross-level) rail inputs have been obtained in a ten- 
degree-of -freedom model. Elastic properties of the passenger car body were repre- 
sented by the bending and torsion of a uniform beam, and the rail-to-car (truck) suspen- 
sion systems were modeled as simple spring-mass -dashpot oscillators. Lateral spring 
nonlinearities approximating certain complicated truck mechanisms were introduced and 
their effects on car response were explored. The models have been excited by displace- 
ment and, in some cases, velocity inputs from the rails, and a wide variety of deter- 
ministic (including sinusoidal) and random input functions have been used. Results were 
obtained in the form of acceleration responses in various parts of the models, particular 
attention being given to car responses. From this study, the following observations are 
made. 

A trial-and-error optimization procedure for minimizing acceleration responses 
in the car has been shown to be feasible for a dynamic system with a few spring and 
damping elements. This feasibility was demonstrated for the vertical model, in which a 
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systematic variation of stiffness and damping properties showed that larger optimum car 
suspension damping was required for random than for sinusoidal inputs and that it was 
not possible to choose values of system parameters that simultaneously minimized 
accelerations in the car and in other parts of the model. The lateral model contained 
too many spring and damping elements for a practicable optimization by this approach 
and would require instead the application of a computer algorithm. 

The stiffness of the equalizer in the truck suspension was found to be important in 
determining optimum damping in the vertical car suspension (bolster) but was not impor- 
tant for the optimum equalizer damping nor the optimum stiffness and damping of a 
simulated transformer mounting beneath the middle of the car. 

Knowledge of car stiffness and weight properties is needed so that with proper 
design its natural frequencies are sufficiently removed from other frequencies in the 
system to avoid resonances characterized by high car acceleration responses. Torsional 
stiffness properties appeared to be less important than either vertical or lateral bending 
stiffness. Structural damping in the car was shown to be important at car bending and 
torsional frequencies within range of other frequencies in the system. However, car 
responses at car frequencies well removed from this resonant region were insensitive to 
variations in car structural damping. Car responses were also shown to be relatively 
insensitive to wide variations in car weight. 

Car acceleration responses to sinusoidal input (transfer functions) were shown to be 
sensitive to variations in vertical and lateral bolster damping and relatively insensitive to 
variations in transformer stiffness (frequency) and damping. 

Rolling stiffness and displacement of the traction motor in the truck suspension 
were shown to have minor effects on car acceleration responses in the lateral model. 

Results of the study for the lateral model indicate that the real-time simulation 
(RTS) method could be used in preliminary design and troubleshooting of rail car config- 
urations. Wide differences in lateral car response were shown for different car speeds 
and the various transient and random inputs that were available with this method, and 
results further showed that no fixed value of lateral bolster damping would reduce lateral 
car accelerations at all car speeds or in all car locations. 

Lateral car acceleration responses were higher for nonlinear than for linear springs 
for two different inputs, and these responses are shown to be particularly sensitive to 
variations in lateral bolster damping. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., June 30, 1971. 
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APPENDIX A 


DERIVATION OF EQUATIONS OF MOTION FOR VERTICAL MODEL 


This appendix presents the derivation of the matrix elements of equations (1) and (4) 
for the vertical model. These equations were obtained from Lagrange's equation in the 
form 


d/6T\ 3D 3U 
dt \3qj] + + Sq j 


0 


(Al) 


where T is the kinetic energy which is assumed here to be a homogeneous quadratic 
function of the velocities, D is the dissipation energy, U is the potential energy, 
and qj is a generalized coordinate. These energies are given in terms of the degrees 
of freedom of the mathematical model of figure 1. Pitching and rolling motions of the car 
are neglected, and small displacements are assumed throughout the model, consistent 
with linear theory. Performance of the indicated operations in equation (Al) for each 
degree of freedom leads to a set of simultaneous equations of motion represented by the 
general matrix equation (4). 


Energies of the Vertical Model 

Kinetic energy .- The kinetic energy of the four-degree-of-freedom model in figure 1 
is given by 

T v = y ffot) + z c ^3 2 dx + f M g (f g( l 3 2 + m A |zA(tjj 2 (A2) 

where m is the mass of the car per unit length L, Mg is the mass of the transformer, 
and M^ is the mass of the traction motor in each truck. The bending deformation of the 
car is given by 

w(x,t) = a(t) W(x) (A3) 

where a(t) is the generalized coordinate having the units of displacement and W(x) is 
the mode shape. The vertical traction-motor displacement z^ is the same for both for- 
ward and rear trucks. 
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APPENDIX A - Continued 


Pot ential energy .- The potential energy of the vertical model is 



dx + 


-jjz(d,t) - z A (t| 2 + jz(L - d,t) - z A (t]| 2 J 


+ Y - z g(§ 2 + k A Ea 1 W “ S(§ 2 (A4) 

where (EI) y is the bending stiffness of the car in the vertical direction, z(d,t) is the 
vertical displacement of the rear trucks, z(L - d,t) the vertical displacement of the 
forward trucks, z^,tj * s the vertical displacement of the middle of the car, and S(t) is 
the displacement input from the rail to both trucks. The truck and mid-car displacements 
are related to the four degrees of freedom by 


z(d,t) = w(d,t) + z c (t) = a(t) W(d) + z c (t) 

z(L - d,t) = w(L - d,t) + z c (t) = a(t) W(L - d) + z c (t) 

z d’ 1 ) - "(H + z c<*> = a <‘> w (f) + z c<‘> 
Substitution of equations (A5) into equation (A4) gives 


•"V 



u v = 


(El), 


|a(t) W(d) + z c (t) - z A (§ 


+ |a(t) W(L - d) + z c (t) - z A W] 2 j 
+ ^[a(t) + z c (t) - z g (t>] 2 + k A jz A (t) - S(§ 2 


(A5) 


(A6) 


Diss ipation energy .- The dissipation energy is given in terms of viscous and struc- 
tural damping constants by 


Cc 

D - — 
u y 9 


|ji(t) W(d) + z c (t) - z A (tj ) 2 + JI(t) W(L - d) + z c (t) - 



+ 





2 


+ 


o GO ( EI )v 



dx 


(A7) 
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APPENDIX A - Continued 


where g cy is the structural damping coefficient of the car in vertical bending. By 
following the approach of reference 8, structural damping is introduced in equation (A7) 
with its magnitude proportional to the elastic restoring forces in the car and in phase with 
the velocities of oscillation of the car. The presence of the frequency co in the denomi- 
nator of the last term of equation (A7) limits the use of the structural damping concept to 
systems with sinusoidal motions. In this investigation, structural damping was included 
only in the acceleration transfer- function computer program for equations (4). 

Car Bending Mode 

The bending mode shape of the car in equation (A3) is taken as the bending mode of 
a uniform free- free beam and is given by 

W(x) = W m (x) = cosh /3 m x + cos j3 m x - Q' m (sinh j3 m x + sin |3 m x) (A8) 

where j3 m and Q! m are obtained from the equations 

cosh /3 m L cos /3 m L = 1 

cosh /3 m L - cos 0 m L 

Q^m = * — 

sinh j3 m L - sin /3 m L 

for m = 1, 2, . . . Solutions of equations (A9) are given in various places in the liter- 
ature. (See, for example, ref. 9.) Since it was assumed that the motion of the car could 
be adequately represented by using only the fundamental mode, it was sufficient to con- 
sider m = 1, for which /3jL = 4.730 and aj = 0.9825. This assumption is essentially 
the same as that made in reference 4, except for the different approximation made for the 
mode shape. Mode- shape values for different stations along the car are given in table VII 

Matrix Elements of the Mathematical Vertical Model 

The elements of the matrices [m], jjc], and [k] of equations (4) were obtained 
by carrying out the operations indicated in Lagrange’s equation (Al) for each degree of 
freedom in the model and these elements are listed in tables VIII and IX. All matrix ele- 
ments of equations (4) were obtained by dividing through each of the equations by the car 
mass M c and the sinusoidal input amplitude S Q . The elements of ^q-|J>, 
and (Qo\ are as follows: 
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In the matrix elements, the natural frequencies of various parts of the model con- 
sidered separately (a>B> w s» w g> and the dimensionless viscous damping coef- 

ficients, are related to the spring mass and damping constants by 
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First vertical bending frequency of car: 


= 03lL) 2 ' (EI)v 


1 M C L K 


(Alla) 


Bolster rigid-body vertical translation frequency: 



Transformer frequency: 



(A lib) 


(Allc) 


Equalizer spring- traction- motor frequency: 
co A = 

Bolster damping coefficient (for both trucks): 

- c s 
^ " M c o> s 

Damping coefficient for the transformer suspension: 

- c g 

” 2MgWg 

Damping coefficient of equalizer suspension in each truck: 



C A 

2M A co A 


(A lid) 


(Alle) 


(Allf) 


(Allg) 


where M c = mL is the mass of the rail car. The equation for cog is a reduced form 
of the equation 



(A12) 
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APPENDIX A - Concluded 


Using equation (A8) for W m (x) leads to W m 2 (x) dx = L and 

r L [d 2 W m (x)] 


dx 2 


dx = L * Substitution into equation (A12) of these values of the integrals 


leads to equation (Alla). It should also be noted that in the derivation of equations (4), 


f L 

there is an off-diagonal element Mj 2 = m \ W m (x) dx. However, by virtue of conserva- 

J 0 

tion of momentum in a free-free beam, this integral vanishes, and the mass matrix 
reduces to a diagonal matrix. Evaluation of the three foregoing integrals is found in 
reference 10. 
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APPENDIX B 


DERIVATION OF EQUATIONS OF MOTION FOR LATERAL MODEL 

This appendix presents the derivation of equations (1) and (4) for the lateral model. 
These equations were obtained in the same manner as those for the vertical model, 
beginning with Lagrange's equation (eq. (Al)). The derivation of equation (1) also involved 
the introduction of terms into the potential energy to allow for the presence of nonlinear 
lateral bolster springs, and these terms were different from those corresponding to the 
linear lateral bolster springs. The lateral model is further assumed to be uncoupled 
from the vertical model. 


Energies of Lateral Model 

Kinetic ener gy.- The kinetic energy of the ten-degree-of -freedom model shown in 
figure 13 is: 


= y J [v(x,t) + y c (t) - h 4 0 c (x,t) - _ x \ ?>(t) dx 



dx + | M g [y g (t )] 2 


h M A 



2 


+ i M Af [>Af ( fc )] 2 + hr pAr (t) ] ' + Kf [*Af wf < B1 ) 

where I c is the mass moment of inertia of the car cross section in roll about the section 
center of gravity. The lateral bending deformation of the car is given by 

v (x , t) = b (t) W (x) (B2) 

with b(t) the generalized coordinate for the car and W(x) the bending mode shape 
given by equations (A8) and (A9). Rigid -body roll and torsion of the car are contained in 
the following expression for 6> c (x,t): 

e c (x ,t) = e c0 (t) + e c t (t) cos ^ (B3) 
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APPENDIX B - Continued 

Potential energy .- The general expression for the potential energy of the lateral 
model with linear springs may be written as 


U L= — 


(ei) t . r L a 2 


GJ f L ®0 c (s,t) J 

dx + — - \ — dx 

2 Jq 9x 


+ ~ 1 t \ y A r (t) - ~ h 2rV (t) “ h 3%(d,t) 


+ -f^ljAf^ - y( L - d >t) - h 2f e Af (t) ' h 3 0 c< L - M 

k r / \ / v 1 2 

+ _ ?[_ y g (t) ' y (l’7 " h T e ci^> t ) 


+ t/[| ' 9c(d ’ t) " I e Ar (t) l + F I + f 


+ |e c (L - d,t) - f ^(t) + - |a c (L - d,t) + f ^(t) 


+ f V (t) - S r(t) - 6 r(t) + - f e Ar (t) - S r (t) 


+ | 0 Af (t) - S f(t) - 6 f(t) + - 1^(0 - S f (t) 


T [ y ArW + h lr e Ar<# ' Y r®] 2+ ^ [^Af® + h lf 0 Af^ ‘ Y f®] 


k «r a 2 ... k <vf a 2 ... 
+ — 0 Ar<*> + — 0 Af ® 


37 



APPENDIX B - Continued 


The lateral and torsional displacements of the car at the trucks and the middle of the car 
are given by 


y(d.l) = v(d,t) + y c (t) - ^ - dj 77(1) = b(t) W(d) + y c (t) - (y - d) n(t) 

y(L - d,t) v(L - d,t) + y c (t) + (y - d ) >7(0 = b(t) W(L - d) + y c (t) + - dj ??(t) 

y (l' 1 ) ■ V (H + y c ( ° = b(t) W (l) + y c« 

0 c (d,U - ti c0 (t) 4 » cl (t) cos Hi 

«c< L “ d.t) - " c q ( t ) 4 (t) cos — — — — f 


(B5) 


Substitution of equations (B2) , (B3), and (B5) into equation (B4) gives the following 
expression for potential energy as a function of the generalized coordinates: 


u (ei) l r L 

u l=-t-3 0 


b(t) 


d 2 Wfc) 


dx z 


dx + 


GJ f ^ 7 x 7 T . 77 X 

— J 0 [_ cl( ) L sin r 


dx 


H cr 


y Ar (0 - b(t) W(d) - y c (t) 4 . dj tj(t) . h 2r e Ar (t) - li 3 ^d c0 (t) + fi cl (t) cos 


+ "ir^Af^ - b(t) W(L ' d) - y c (t) - (| - d ) bft) - h 2I e M (t) - h 3 


y c0® + fi cl® cos 


!T(L - d) 




y g (t) - b(t) - y c (t) - h T 0 cO (t) 


fr 


e c ()W + e clW cos T ' e Ar (t) 


0 cO (t)4e cl (O cosiO^d) -e M (t) 


-,2 


k Ar J [~R 

4 L 2 


e Ar (t) - s r (t) - e r (t) 4 


|e Ar (t) + s r (t) 


Af 


-,2 


|d At (t) - S f (t) - 6 t (t) 


-.2 


4 S f (t) 


Iv 


-.2 


y Ar (‘) + h lr e Ar (t) - Y r (t) 


y Af (t) 4 - Y f (t) 


k £»r ,2 , > kjyf 0 , . 
— e Ar W + — 9 AfW 


(B6) 
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APPENDIX B - Continued 

With nonlinearities allowed in the lateral bolster springs, kj cr and are no 

longer constant, and the potential energy terms associated with these variable spring con- 
stants must therefore be replaced. The replacement terms are derived by first defining 
a variable spring function k(£) such that the restoring force is k(£) £, where £ is 
the spring deformation. The work done by this force moving through a differential 
displacement d£ is [k(£) IJdi;, and the potential energy developed in the spring 
stretched to some displacement Ay is then obtained from the integral 

%L - |j(« 1 « f 67 ' 

The displacement Ay is a function of the degrees of freedom of the system, and in this 
investigation, Ay is given by equations (15). Thus, the potential energy due to the non- 
linear lateral bolster springs is 

%l - J 0 iyr f ter«> g + J 0 ay£ gfcf «> g« <“> 

which replaces the second and third terms in equation (B6). 

Dissipati on energy . - The dissipation energy of the lateral model in terms of both 
viscous and structural damping constants is written as follows: 


) - te 


y Ar - b(t) W(d) - y c (t) + - dj 7/(t) - h 2r f? Ar (0 - h 3 ^ c0 (t) + 0 cl (t) cos ^ 


y M (t) - b(t) W(L - d) - y c (t) - ^ - dj ^(t) - h 2f 0 Af (t) - h 3 ^ c0 (t) + 0 cl (t) cos 

2 


tt(L - d) 
L 


C j£ 

2 

CoB 2 


u Ar 


Af 


y g (0 - b(t) w(tj. y c (t) - h T e c0 (t) 


e cO (t) + 0 c icos^-e A r(t) 


e c0 (t) + e cl cos 7r(L . d) - e Af (t) 


—i 2 


| a Ar (t> - s r (t) - o r (t) 


->2 


| e Af (t) - Sf(t) - s f (t) 


| e Ar (t) + s r (t) 


f e Af (t) + s f (t) 


~lr 




}*,«) * »lr»Ar<‘> - Y r«> 


-i 2 


y M (t) + h lf 0 Af (t) - Y f (t) 


b(t) 


d^W(x) 


—i 2 


dx^ 


dx 


i g cT rT r 1 
2— GJ 1 


-.2 


.6 cl (t) l sin f 


dx 


(B9) 
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where g c j is the structural damping coefficient of the car in lateral bending, and g c rp 
is the structural damping coefficient of the car in torsion. These coefficients have been 
introduced in the same manner as the bending structural damping coefficient of the car in 
the vertical model. (See appendix A.) 


Matrix Elements of Mathematical Lateral Model 

The elements of the square matrices of equations (1) and (4) were obtained by 

carrying out the indicated operations in Lagrange's equation for each degree of freedom 

in the lateral model in the same manner as for the vertical model in appendix A. For 

the nonlinear springs, the Lagrange operation reduces the integrals in equa- 

8q. 

tion (B7) to functions of Ay; that is, 3 


8U, 


NL 


8m 




(BIO) 


Thus, the integrals in equation (B8) reduce to 


8U, 


NL 

9 <h 


= [vM K) + [>f H] M 


(Bll) 


where qj = b(t),y c (t), . . . 0^(t). The elements of the stiffness matrix obtained from 
this operation for equation (1) are identical in form with those obtained for the linear 
springs on the basis of equation (B6). 


Matrix elements for the lateral model are given in tables X to XII. Each matrix is 
symmetric about its principal diagonal. The mass matrices are given in table X, the 
upper matrix being for the general form of equation (1) and the lower matrix, for the 
frequency-domain form of equations (4). The stiffness and damping matrices are given 
in table XI for equation (1) and in table XII for equations (4), the upper matrices being the 
stiffness (K) matrices and the lower matrices the damping (C) matrices. The moment 
of inertia of the car per unit length about the car elastic axis is given by I^a = Ic + nih|. 
The matrix elements corresponding to equations (4) were obtained by dividing each equation 
by the car mass M c (or mL) and sinusoidal input amplitude Q 0 in a manner similar 
to that done for the vertical model in appendix A. 


and k 


Ic f 


in 


The rearward and forward lateral bolster spring constants (k^ cr 
table Xl) were kept equal throughout the investigation. However, the distinction between 
them is necessary when these springs are nonlinear, in order to relate the nonlinear 


40 



APPENDIX B - Continued 


deflections properly to the motion throughout the system in accordance with equations (15) 
during the numerical integration of equation (1). Since there is no nonlinear damping 
in the system, c^ cr = c^ = c^, in the damping matrix of table XI. 

The input matrices of equation (1) for the lateral model are 


{ 8 } ■ 


'0 


r 0 

1 

1 . 

* 

0 

II 

° 1 

k, r Y r (t) 


c * r Y r (t) 

k K Y f (t) 


c Zf Y f (t) 

k Ar R 

k ir h lr Y r(t) + 6 r« 


c Zr h lr Y r« + 6 r^ 

k »rR 

[kj f hi f Y f (t) + -^-6 f ( t ) J 


c ff h lfYf (t) + ^ 6 f (t) 

v y 


(B12) 


Equations (4) were programed for either lateral or cross -level rail inputs, and the input 
amplitudes Y q j and 6 0 f to the forward trucks were chosen as reference displacements 
in defining the dimensionless degrees of freedom for the {Oq} and matrices. 

The complete sinusoidal relations involving these input amplitudes are 


Y f (t) = Y of sin cot "1 


8 { (t) = 6 o{ sin wt I 

Y r (t) =Y of sin (wt + KjU 
Y f (t) = 5 0f sin (cot + k 2 'J 


(B13) 


where and k 2 are phase angles which, in this study, were each an integral multiple 
of it. Even multiples imply in-phase inputs to front and rear trucks, whereas odd 
multiples imply out-of-phase inputs. The quantity S r (t) or 6j(t) represents the 
difference in vertical inputs to left and right rails under each set of trucks (see again 
fig. 13(b)) and completely defines the cross-level input. The quantities S r (t) and Sf(t) 
disappear in the derivation of the equations of motion. The matrix elements for , 

{^2}’ (Qi}> and {^2} for these inputs are as follows: 
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' b l " 


b 2 

Y of 


Y of 

y c i 

"■ 

y c2 

Y of 


Y of 

e c01 B 


e c02 B 

Y of 

I 

Y of 

e cll B 


e cl2 B 

Y of 


Y of 

J7lL 

r . 

V2 h 

< Y of 

> {* 2 J = \ 

Y of ) 

Zsi 


M 

Y of 


Y of 

^Arl 


y Ar2 

Y of 


Y of 

e Afl R 


e Af2 R 

Y of 

k J 


Y of 

J 


0 


0 


Q 


'l = co? 

J lr M,, 


2 M Ar sin (cot + kQ 


sin cot 


co 


2 M Af 
if ML 


CO 


M Ar h lr sin(cot+K 1 ) 


ir M c R 


sin cot 


CO 


2 M Af h lf 
if M c R 


q 2 = 


2n 7 w 




M Ar cos ( wt + k 2) 


cos cot 


M Af 

2n 7f co 7f 

U li m c 


2n (jo ^Ar h l r C ° S ( C ° t + ^ 2 ) 
i T i r Me R cos cot 


hif 

2 ”«“« 


J 


) (B14) 
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for the lateral rail input. The {Qj} and {^2} matrices for the cross-level input 
are of the same form as those for the lateral input except that 6^ replaces Y^. The 
cross-level input matrices are 




0 1 

CO 2 lAr 

sin (cot + 

Ar M C R 2 

sin cot 

w Af 

*Af 

MgR 2 


Q2H 


> (B15) 


: Ar cos(ft>t+ k 2 ) 


Ar Ar McR 2 cos 


2n Af w Af 


*Af 

McR 2 




As in the case of the vertical model, the natural frequencies and viscous damping 
coefficients of the lateral model appearing throughout the matrix elements in tables X 
to XII and equations (B14) and (B15) are related to the spring and damping constants by 
the following expressions: 


First lateral car bending frequency: 


U) 


IB 


“ (% L ) : 


2 FD] 


MfL' 


Lateral shear bolster rigid -body translation frequency: 


^c = 


k fc r + k Jcf 


M r 


Bolster rigid-body rolling frequency: 
= B. 


O) 


oT 


*-s 


2I EA l 


First torsional frequency of car: 


_ _GJ_ 
T L 'IlEA 


, , _ 7T 

( a ) m — 


(B16a) 


(B16b) 


(B16c) 


(B16d) 
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Lateral transformer frequency: 


0), = I 

lg \j 


be. 


M, 


g 


Lateral equalizer frequency: 

V 




m a 


Rear equalizer rolling frequency: 


“Ar = f 


Ar 
\|M Ar 


Front equalizer rolling frequency: 


W Af “ f 


Af 


PAf 

Roll frequency of rear traction motor: 


co 


ar 


ar 


\|I A r 

Roll frequency of front traction motor: 


U) t = 

a f 


J»f 

[ Af 


Bolster lateral damping coefficient (for both sets of trucks): 
c, 


n,„ = 


te 


Lateral transformer damping coefficient: 
C i 


n ^g 2w^gMg 

Bolster rigid-body roll damping coefficient: 


1&. 


CoB^ 


n 


oT 4 o> oT I ea L 


(B16e) 


(B16f) 


(B16g) 


(B16h) 


(B16i) 


(B16j) 


(Bl6k) 


(B161) 


(B16m) 
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Lateral shear equalizer damping coefficient: 
n - "A 

l 2o^M a 

Rear equalizer roll damping coefficient: 


n. 


c Ar R2 


Ar 

Front equalizer roll damping coefficient: 
c Af R2 


” M = 

Roll damping coefficient of rear traction motor: 


n = 




«r 2 co ffr I Ar 
Roll damping coefficient of front traction motor: 
c oti 


= 


ai ZUaftM 


(B16n) 


(B16o) 


(B16p) 


(B16q) 


(B16r) 
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TABLE I.- PROPERTIES AND INPUTS FOR VERTICAL MODEL 


(a) Fixed properties 


Property Value 

Car length, L, m (in. or ft) 25.9 (1020 or 85) 

Truck attachment distance, d, m (in. or ft) 3.81 (150 or 12.5) 

Transformer weight, M g g, N (lb) 57 800 (13 000) 


Traction motor weight, M^g, N (lb) 53 000 (11 908) 

Bolster (car) rigid-body frequency, f s , Hz 1.0 


(b) Variable properties 


Figure 

! 

Car weight, M c g 

Car bending 
frequency, 
f B , Hz 



Transformer 
frequency, f„, 
Hz g 

Equalizer 
frequency, 
f A , Hz 

Car damping 
coefficient, 
gcv 

Bolster 

damping 

coefficient, 

n s 

Transformer 

damping 

coefficient, 

n g 

Equalizer 

damping 

coefficient, 

% 

Input 

kN 

lb 

3 

574 

129 000 

8.0 

6.0 

4.12 

0 

, 

0.11 

0.25 

0.10 

Sinusoidal 

j 4 

574 

129 000 

8.0 

Variable 

4.12 

0 

.11 

.18 

.10 

Sinusoidal 

5 

574 

129 000 

8.0 

6.0 

4.12 

0 

.11 

Variable 

.10 

Sinusoidal 

6 

574 

129 000 

8.0 

6.0 

4.12 

0 

Variable 

.25 

.10 

Sinusoidal 

7 

574 

129 000 

8.0 

6.0 

4.12 

0 

.11 

.25 

Variable 

Sinusoidal 

8 

574 

129 000 

8.0 

6.0 

4.12 

0 

Variable 

.25 

.10 

Random* 

9 

574 

129 000 

8.0 

6.0 

4.12 

0 

.30 

.25 

Variable 

Random* 

10 

574 

129000 

8.0 

6.0 

4.12 

0 

.30 

.25 

.30 

Random* 

11 

574 

129 000 

Variable . 

5.0 

4.12 

Variable i 

.057 

.18 

.03 . 

Sinusoidal 

12 

Variable 

Variable 

Variable 

| 5.0 1 

4.12 ■ 

i 0 ! 

.057 

1 - 18 

I - 03 

1 Sinusoidal 


* Spectral density given by equation (3). 



TABLE n.- OPTIMUM PARAMETER VALUES 
FOR TWO EQUALIZER FREQUENCIES 


Parameter 

Values for equalizer 
frequency, f A , of — 

4.12 Hz 

5.53 Hz 

Transformer frequency, f g , Hz . . 

6 (or 6.5) 

6 

Transformer damping, n g 

0.25 

0.18 

Bolster damping, n s 

0.11 

0.057 

Equalizer damping, n A 

0.10 

0.10 


TABLE HI.- TRANSFER FUNCTIONS FOR LATERAL MODEL 


Location 


Car 


Transformer 


Traction 

motors 


In car over 
bolsters 
(vertical) 


Equalizers 

(vertical) 


Acceleration transfer function 


A cl (x) 

S% 


b^Oc) + y cl - y(®c01 + %11 cos S) 





b 2 W(x) + y c2 - 



+ &cl2 cos 


7 TX\ 

L ) 



Kl w 2 fl2 IF 
^ = -fgl + y g2 


A Af _ w 2 I -2 -2 

i Q^- T ^Afl + y A f2 


t*L = /yf 1 + y 2 

gQo g r Arl Ar 


A Bf _ 

gQ 0 2g 


n Br _ 

gQ 0 ~ 2g 


%01 + %11 cos 


^c02 + %12 cos Tjj T rd) 


2 2 
(%01 + ®cll cos lj) + (%02 + ®c!2 cos L~) 


L 

1/2 


1/2 


Ef _ w* s2 , b2 


g Q: = 2FrAfl + e Af2 


A Er 

gQ 0 




fl 2 4. fl 2 

e Arl + y Ar2 




TABLE IV.- FIXED PROPERTIES FOR LATERAL MODEL 


Property 

Car length, L, m (in. or ft) 

Truck attachment distance, d, m (in. or ft) 

Lateral distance between bolster springs, B, m (in.) 

Lateral distance between equalizer springs, R, m (in.) 

Distance of traction motor to lateral equalizer, h-^ r and h^, m (in.) 

Distance of traction motor to lateral bolster, h 2 r and h 2 f, m (in.) 

Distance of lateral bolster to car elastic axis, h 3 , m (in.) 

Distance of lateral transformer to car elastic axis, h^p, m (in.) 

Distance of car elastic axis to center of gravity, h^, m (in.) 

Center-to-center distance between truck axles, (AZ) t , m (in.) 

Transformer weight, Mgg, kN (lb) 

Traction-motor weight, M Ar g and M^g, kN (lb) 

Traction-motor mass moment of inertia, I Ar and Ipj, kg-m^ (in. -lb-sec . . . . 

Lateral bolster spring constant, k^j. and k^ cf , N/cm (lb/in.) 

Vertical bolster spring constant, k s , N/cm (lb/in.) 

Lateral equalizer spring constant, k^ r and k#, N/cm (lb/in.) 

Vertical equalizer spring constant*, k Ar and k^, N/cm (lb/in.) 

Lateral bolster frequency, f^, Hz 

Bolster rolling frequency, f 0 rp, Hz 

Lateral equalizer frequency, f^ r and f^f, Hz 

Equalizer rolling frequency, f Ar and f^, Hz 

Equalizer roll damping coefficient, n Ar and n ^ 

Traction-motor roll damping coefficient, = n^f 

* Total of four springs per truck at 9070 N/cm (5 170 lb/in.) each. 


Value 


25.9 (1020 or 85) 

3.81 (150' 

or 12.5) 

2. 

34 (92) 

2. 

01 (79) 

0.534 (21) 

0.305 (12) 

0.305 (12) 

0.305 (12) 

0.610 (24) 

2.595 (102) 

57.8 

(13 000) 

53.0 

(11 908) 

2100 

(18 600) 

5 250 

(3 000) 

11520 

(6 570) 

7 000 

(4 000) 

36 300 

(20 680) 


0.674 


0.997 


1.81 


6.62 


0.10 


0.60 


TABLE V.- VARIABLE PROPERTIES 


Figure 

Traction-motor 
rolling frequency, 
for ~ ^orf’ 

Hz 

Car lateral bending 
frequency, f iq , 

Car torsional 
frequency, f-p, 

Lateral transformer 
frequency, fjg, 

Car weight, 

Car mass moment 
of inertia per 
unit length, IgA 

Car structural damping 
coefficient in lateral 
bending, g ci 

Car structural 
damping coefficient 
in torsion, g c p 





kN 

lb 

kg-m 

lb-sec^ 

19 

18.0 

8.0 

15.0 

6.0 

574 

129 000 

3090 

695 

0 

0 

20, 21 

Variable 

8.0 

15.0 

6.0 

574 

129000 

3090 

695 

0 

0 

22 

18.0 

Variable 

15.0 

6.0 

574 

129 000 

3090 

695 

Variable 

0 

23 

18.0 

8.0 

Variable 

6.0 

574 

129000 

3090 

695 

0 

Variable 

24 

18.0 

Variable 

Variable 

6.0 

Variable 

Variable 

Variable 

Variable 

0 

0 

25 

18.0 

8.0 

15.0 

Variable 

574 

129 000 

3090 

695 

0 

0 

26 

18.0 

8.0 

15.0 

6.0 

574 

129000 

3090 

695 

0 

0 

27 

18.0 

8.0 

15.0 

6.0 

574 

129 000 

3090 

695 

Variable 

Variable 

28 

18.0 ! 

8.0 

15.0 

6.0 ! 

574 

129000 

3090 

695 

0 

0 

29 

18.0 

8.0 

15.0 j 

6.0 

574 

129000 

3090 

695 

0 

Variable 

30 ] 

18.0 ! 

8.0 

15.0 

6.0 

574 

129 000 

3090 

695 

0 

0 

1 

31, 36 ! 

18.0 , 

1 

8.0 

| 

15.0 

6.0 ! 

574 

129000 

3090 

695 

0 

0 

32 

18.0 

8.0 | 

15.0 

6.0 - 

574 

129 000 

3090 

695 

0 

0 i 

33 j 

18.0 

8.0 

15.0 | 

6.0 

574 

129 000 | 

3090 

695 

0 

o t 

34 

18.0 

8.0 j 

15.0 | 

6.0 

574 

129 000 | 

3090 

695 ; 

o ; 

o 

35 

18.0 


15.0 

6.0 

574 ! 

129 000 [ 

3090 j 

695 

0 i 

0 | 

37 , 

18.0 

8.0 j 

15.0 

1 

6.0 l 

574 i 

129 000 ! 

3090 j 

695 

o i 

0 

38 ! 

18.0 

8.0 

15.0 

6.0 1 

574 | 

129 000 j 

3090 ! 

695 ; 

° i 

o ! 

39 ' 

18.0 ! 

8.0 

15.0 

6.0 

574 

129000 1 

3090 : 

695 

o 

0 


AND INPUTS FOR LATERAL MODEL 


Lateral transformer 
clamping coefficient, 

Lateral bolster 
damping coefficient, 
n Zc 

Bolster roll 
clamping coefficient, 
n oT 

Lateral equalizei 
damping coefficient, 
n Zr = n Zf 

Sinusoidal lateral 
input (fig. 15(a)) in-phase, 
out-of- phase: Kj = 0 or it 

Sinusoidal cross -level 
input (fig. 15(a)) in-phase, 
out -of -phase: k^~ 0 or tt 

Transient and 
random inputs 
(figs. 15(b) to 15(g)) 

0.25 

0.057 

0.20 

0.10 

0 



.25 

.057 

.20 

.10 

0, -- 

n 


.25 

.057 

.20 

.10 

0 

— - 


.25 

.057 

.20 

.10 

■n 

TT 


.25 

.057 

.20 

.10 

JT 

J T 


.25 

.057 

.20 

.10 

0 



Variable 

.057 

.20 

.10 

0 

___ 


.25 

Variable 

.20 

.10 

0, JT 

___ 



.25 

1 

.20 

.20 

Variable 

0, IT 

— 


.25 1 

1 .057 

Variable 

..0 

' 

0, JT 


.25 

Variable 

Variable 

.10 

— 

— 

fig- 150)) 

.25 

Variable 

.20 

.10 

— 

___ 

fig- 15(c) 

.25 

Variable 

.20 

.10 

--- 


fig- 15(d) 

.25 

.057 

.20 

.10 

— 

1 

fig. 15(d) | 

.25 

.057 

.20 

.10 


| 

fig. 15(c). fig. 15(f) 

.25 

Variable 

.20 

.10 



fig- 15(g) 

.25 

Variable 

.20 

1 

.10 


___ 

fig. 15(e) 

.25 

Variable 

.20 

.10 


... 

fig. 15(e) : 

.25 

Variable 

.20 

.10 


... 

fig- 15(d) 
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Ill I III I 


I I II I 


TABLE VI.- ROOT-MEAN-SQUARE ACCELERATIONS 
FOR TWO DIFFERENT LATERAL RANDOM INPUTS 
TO LATERAL MODEL 



RMS accelerations, g, for - 

Car location 

No. 3 

Eastern U.S. railroad 
test track (fig. 37) 

Combined input 
(fig. 35) 

Rear trucks 

0.0416 

0.0342 

Forward trucks 

.0416 

.0342 

Middle of car 

.0465 

.0334 


TABLE VTL- CAR BENDING MODE SHAPE 


X 

L 

w i(e) 

0 

2.000 

.1 

1.0743 

.147 

.6509 

.2 

.1954 

.3 

-.5440 

.4 

-1.040 

.5 

-1.2156 

.6 

-1.040 

.7 

-.5440 

.8 

.1954 

.853 

.6509 

.9 

1.0743 

1.0 

2.000 
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TABLE VHI.- MASS MATRIX FOR VERTICAL MODEL 


1 2 

3 

4 

1 1 0 

0 

0 

2 0 1 

0 

0 

3 0 0 

Mg 

M c 

0 

4 0 0 

0 

2 ^A 

M c 
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TABLE IX.- 


,2r 


W z (d) + W^{L - d) 


+ ai 2 ^£-W 2 (t 

g Mc \ 2 


jFrom table VII 
2 


-j~W(d)+W(L-d)j 

A 


gMc 


Wit 




ol-& 


g Mc 


n s“sR 


W 2 (d)+W 2 (L-d) 


+ 2 ^TT^ 2 k) + Sr 


Me 


3 CV C 0 



r "| 


n s“s 

W(d) +W(L -d) 

2 ( ns " s+n g^'“g) 

+ 2n g 

Ss w (i) 


Mg 

- 2n g^“g 


-n 0 w = rw(d) + W(L-d)l 


- 2110(00 


MODEL 


L2_^ 

, W(d) = W(L - d>] 



0 


2 ln s (o s + 2n 





m 




TABLE XI.- GENERAL STIFFNESS AND DAMPING 


(El) L /3jL + k ter W 2 (d) 

+ k fcf W 2 (L-d)+k ig w2^j 


k fcr W(d) + k fcf W(L-d) 

+k fe w (r) 


h 3[ k fcr W(d)+k tef W(L " d) ] 

+ Vt w (?) 


kfcrWWcosSl 


+ k tef W(L -d) cos ~ 


k fcr + k Zcf + k fe h 3( k ter + k tef) + k Zg h T h 3^ k 


ho k,„^ cos + k, r cos - 


3 K Zcr cos "l" + K Zcf 


(2 - d ) [? k fc rW<<,) * k fcf w<L • 


2 ' d i \' k lc t + k fcf } 


“V" +h l( k icr +k tef) 


c ]r W z (d)+W z <L-d) 


C|c[w(d) + W(L-dj] 

+c ie w (r) 


l^[co S 2 ^ + cos2ia^] 

] + h 3 [ k fcr cos2 ff 

[ , „_2 7r(L-d)l . GJL (n\ 2 

I +k tef cos L J 2 W 


h 3 (ff ‘ d )(“ k fcr + k fcf) 


h 3\ 2~ ” / j" k fcr COS "l* 


h 3 c Ic[ w ( d ) + w ( L - d ]] 
+ 1 >T^ W (?) 


2h 3 c Zt: +h T c ig 


C.-B 9 2 

~2 __ + 2h 3 c fc + h T c ig 


h 3 c fe[ w < d ) cos r 


h 3Cfc cos T- + cos 




+ cos 2 M 


c,,, #-d| -W{d) + W(L-d) 


h 3 c fc ^- d -cos^ 


8, " c icfW (L - d) 


-h 3 c fc cos 


9 . h 2r c fc W(d> 


C s B , h h 

^1“ ' fc h 2r h 3 y 


' \ 4 ~ C lc l, 2r"3 


c lc h 2i\2 ' d J 


10 h2fCfcW(L-d) 


c s“ K . 1 

"4— ■ c ic h 2f h 3/ 


^•ws)™¥ 
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Mt r 


MATRICES FOR LATERAL MODEL (EQ. (1)), 




TABLE XU - STIFFNESS AND DAMPING MATRICES FOR LATERAL 


1 

2 

3 

^ B + ^[w 2( d)+ w2 {L-d)] 

^[w(d} + W(L-d)] 

^[w(d) + W{L-d3 

‘sfltf 

3 

+ 

1 

+ 45 w ® 

< 45 ^ w ® 



4^-455 




n fc ^[w 2 ( d ) + W 2 (L-d3 


- 4(^5 

+ 2n Is“2gMf w2 (?) +6 ci"^ 



n te“'te[ W<d)+W(L ' d ] 

2 (''ic“-'fc +n te“is 5 ) 


n fc u> fc ^[w<d) + W(L-d^] 

- 2 v^ 55 w (?) 

/ h 3 Mg h T ^ 

2 (Vic^ + I 

i 

/ I EA h 2 

2 ^VVr“^ 2 + V’ic ^2 

Mg h 2 \ 

+ V*gM^ 

- 

v* £[»««• r 

+ W(L - d) cos j 

! 

W 

y n oT-'oT ^ 4 w p) [ cos r 

«- V 

Vfc TrE WWtW(L - d ^ 

i 

° 

i 

i 

0 

' -v 4 w © 

Mg 

_2 Vfe Me 

l 

' M„ h T 

1 “ 2 Vfe M^ “b 

h 

-"te“ic w < d > 

"W 

_ h 3 

n ic a ic B 

- n te w te W(L-d) 

-Vic 

^*3 

'Vic s' 

w ir w < d > 

n ■ ,l2r 
n ic a Zc R 

/ , %A h 2r h 3^| 

- l n oT u oT m nR " Wc ic* BR ) 

1 .> fc “' te ^W(L-d) 

n w h21 
°ic ic R 

{ , ! EA u h 2f h 3\ 

y n oT u oT m BR " ic ic BR / 


4*3 r, 

2B L 


W(d) cos ^ 


+ W (L - d) cos 




4^r cos ^ +cos ifl^i! 

2 B L L L J 


I 2 *EA ^ 2 h 3 j [" wd 

te B 2 /L COS L 


2\ oT mB 2 


4 )t 


eos^.cos 2 ^ 




I 

'yoT^oT ^g2 


L cos' 


.2 jr(L-d)l 

5 l r 


hg ■ 

Vic b'\ 
s 


+ cos - 


0 


-n 7f> u.’ cos 


ud 

L 




( , ! ea 

\ n oT oT niBR 


, h 2r h 3 l „ nt _ Tri 
- Vic “bIT j C ° S L 


( , ! ea 

~ y n oT a oT niBR 

„ f , h 2f h 3 ) u(L-d) 
' Vic -br7 cos ~ l ^ 


5 



[jw(d) + W(L - d)] 


0 


0 



2 



7rd 

L 


+ cos 


g(L - d)j 



/L 

( 2 


7c icy 


0 



~ n ic“ic 



h /t- 
, h 2r 2 
n lc Ic R y L 


n 2f 

Vic -R- 
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Figure 1.- Mathematical vertical model of railroad car and truck suspension. 










Bolster springs 



Figure 2.- Schematic view of truck suspension system. 





Acceleration transfer function Acceleration transfer function 


CT> 

to 


gA™ 8/in- g/cm g / in . 



Frequency, Hz Frequency, Hz 

(a) End of car. (c) Transformer. 



Frequency, Hz Frequency, Hz 

(b) Middle of ear. (d) Traction motors. 

Figure 3.- Acceleration transfer functions for optimized vertical model. 



I 

I 


g/cm g/in. 


4 — 



(a) End of car. 


3 - 



(b) Middle of car. 


g/cm g/in. 



(c) Transformer. 



(d) Traction motors. 


Figure 4.- Effects of vertical transformer frequency on acceleration transfer functions for vertical model with first 
equalizer spring (f^ = 4.12 Hz). Numbers beside curves are frequencies of peak acceleration transfer functions. 


03 

CO 


05 

g/cm g/in, g/cm g/ in. 



(a) End of car. 


(c) Transformer. 


i.2 r 


w . 0.82 

- ,0^ — o — ^>—0^07.5 


5 . 6 ^^ 


5.0 


O 1 — 


J- 


J 


0 .2 .4 .6 

Transformer damping coefficient, ng 

(b) Middle of car. 



(d) Traction motors. 


Figure 5.- Effects of vertical transformer damping coefficient on acceleration transfer functions for vertical model 
with first equalizer spring (f^ = 4.12 Hz). Numbers beside curves are frequencies of peak acceleration transfer 
functions. 



g/cm g/in. 



Bolster damping coefficient, n s 


(a) End of car. 



g/cm g/in. 



(c) Transformer. 



(b) Middle of car. 


(d) Traction motors. 


Figure 6.- Effects of vertical bolster damping coefficient on acceleration transfer functions for vertical model with 
first equalizer spring (f ^ = 4.12 Hz). Numbers beside curves are frequencies of peak acceleration transfer 
os functions. 

CJ1 


OS 

OS 


g/in. 


g/cm g/in. 


g/cm 



(a) End of car. 




(c) Transformer. 



Equalizer damping coefficient, n^ 


(b) Middle of car. 


(d) Traction motors. 


Figure 7.- Effects of vertical equalizer damping coefficient on acceleration transfer functions for vertical model 
with first equalizer spring (f^ = 4.12 Hz). Numbers beside curves are frequencies of peak acceleration 
transfer functions. 





Figure 8.- Effects of vertical bolster damping coefficient acceleration root-mean-square 
values for vertical model with first equalizer spring ( f ^ = 4.12 Hz). 
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Figure 9-- Effects of vertical equalizer damping coefficient on acceleration root-mean- 
square values for vertical model with first equalizer spring = 4.12 Hz\ 
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Car bending frequency, fg, Hz 

(b) Middle of car. 


Figure 11.- Effect of car vertical bending frequency and structural damping on 
car vertical peak acceleration transfer functions for first equalizer spring 
(f^ = 4.12 Hz). Numbers beside curves are frequencies of peak accelera- 
tion transfer functions. 
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Peak acceleration transfer function 


JL 

in. 


cr 



1 ! J I 

500 550 600 650 700 kN 

Car weight 


Figure 12.- Effects of car weight on vertical peak acceleration transfer function 
for first equalizer spring (f^ = 4.12 Hz) at peak frequencies of about 0.8 Hz. 
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I 


v(x,t) 


X 


Transformer 


Section A - A 


(a) Top and mid-section views. 

Figure 13.- Mathematical lateral model of railroad car and truck suspension, 
with sections A- A, B-B, and C-C viewed from front of car (x = L). 



<3 





s f (t) + 6f(t), S r (t) + 6 r (t) 


(b) Truck suspension (sections B-B, C-C). 
Figure 13.- Concluded. 





-J 


Deterministic 



|(b) Square pulse. 



(c) Square pulse train. 



Random 



Frequency, Hz 

(e) Input spectral densities. 



(f) Random pulse train. 


Q(t,v) 



1 (g) Combined railroad estimate (4) 
and random pulse train. 


Figure 15.- Forward truck inputs. 



Figure 16.- Flow chart for RTS on-line random input capability. DAC denotes digital -to -analog converters; 
ADC denotes analog-to-digital converters; and CRT denotes cathode ray tube. 




Figure 17.- RTS flow chart for solution of equations of motion for 
lateral model. + {o} _ |c]{q} - 








Acceleration transfer function 


g/in. 


g/cm g/in- 



(a) Car. 



(c) Traction motors. 



16 


Frequency, Hz 

(b) Transformer. 



(d) Over bolsters (vertical) 



Frequency, Hz 

(e) Equalizers (vertical). 


Figure 19.- Acceleration transfer functions for lateral model 
with lateral in-phase inputs (/q = 0). 


Acceleration transfer function Acceleration, transfer function 


g/cm g/in. a c1 (O) 

gY 0 f 



0 10 20 30 40 45 


Frequency,. Hz 

(a) End of car. 



Frequency, Hz 


(c) Transformer. 



0 10 20 30 40 45 


Frequency, Hz 

(e) Over bolsters, vertical 
acceleration (rear). 


g/cm g/in. 



0 10 20 30 40 45 

Frequency, Hz 


(b) Middle of car. 



Frequency, Hz 


(d) Traction motor (rear). 



0 10 20 30 4 0 45 

Frequency, Hz 


(f) Equalizers (vertical). 


Figure 20.- Effects of traction-motor rolling degrees of freedom and 
rolling stiffnesses on acceleration transfer functions for lateral 


model with lateral in 


-phase inputs (/q = 0). 
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Figure 21.- Effects of varying traction-motor rolling stiffnesses 
on acceleration transfer functions for lateral model with 
cross-level out-of-phase inputs ^2 = ir )- 
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Figure 22.- Effect of car lateral bending frequency and structural damping on car 
lateral peak acceleration transfer functions. Numbers beside curves are fre- 
quencies of peak acceleration transfer functions. 
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Figure 23.- Effects of car torsional frequency and structural damping on car lateral 
peak acceleration transfer functions for out-of-phase inputs at end of car. Num- 
bers beside curves are frequencies of peak acceleration transfer functions. 
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Figure 24.- Effect of car weight on lateral peak acceleration transfer functions 
at end of car for out-of -phase inputs (ki,/c 2 = ff )* 


83 



i 


g/cm g/in. 




Figure 25.- Variation of car lateral peak acceleration transfer functions with 
lateral transformer frequency. Numbers beside curves are frequencies 
of peak acceleration transfer functions. 
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Figure 26.- Variation of car lateral peak acceleration transfer functions 
with lateral transformer damping coefficient. Numbers beside curves 
are frequencies of peak acceleration transfer functions. 
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Figure 27.- Variation of peak acceleration transfer functions with lateral 
bolster damping. Numbers beside curves are frequencies of peak 
acceleration transfer functions. 
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Figure 28.- Variation of car lateral peak acceleration transfer functions 
with lateral equalizer damping coefficient. Numbers beside curves 
are frequencies of peak acceleration transfer functions. 


86 


Peak acceleration transfer function 


g/cm g/in. 



0 .1 .2 .3 .4 .5 


Bolster roll damping coefficient, n 0 rp 

Figure 29.- Variation of peak acceleration transfer function with bolster 
roll damping for cross -level inputs. Numbers beside curves are 
frequencies of peak acceleration transfer functions. 
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Figure 30.- Lateral car peak acceleration responses to square-pulse 
inputs for X = 0.305 m (12.0 in.). 
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Figure 31.- Lateral car peak acceleration responses to deterministic 

square -pulse train inputs. 
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Figure 33.- Effects of input wave length on lateral car peak acceleration 

responses to wave inputs. 
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Figure 34.- Comparison of lateral car acceleration responses for deterministic 
and random lateral square-pulse train inputs, Y Qr = Y q j = 2.54 cm (1.0 in.), 
deterministic, and root mean square, random. 
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Figure 35.- Lateral responses to combined random pulse train 
and railroad estimate (fig. 15(g)). 
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Figure 36.- Variation of lateral car 
length for deterministic 
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Figure 37.- Lateral car root -mean- square acceleration 
responses to lateral railroad test track input (no. 3 in 
fig. 15(e) with Y or = Y 0 f = 2.54 cm (1.0 in.) rms). 
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Figure 38.- Comparisons of linear with nonlinear lateral car acceleration 
spectral densities for lateral band-limited white-noise input (no. 2 in 
fig. 15(e)). 
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Figure 39.- Nonlinear lateral responses to lateral wave input 
for Y or = Y Q f = 5.08 cm (2.0 in.). X = 9.52 m (375 in.); 
v = 35.8 m/sec (80 mph). 
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